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(57)Abstract: 

PROBLEM TO BE SOLVED: To provide a control 
device for an internal combustion engine for optimally 
controlling the internal combustion engine by utilizing 
the pulsation of pressure in intake pipe for easily . 
estimating a cylinder filled air amount for each 
cylinder. 

SOLUTION: The control device comprises a throttle 
passing air amount calculating means 19 for 
calculating a throttle passing air amount mt for air to 
pass through a throttle valve 18, an extra air amount 
calculating means for calculating an extra air amount vv..^«t^« 
for each cylinder equivalent to a drop amount A 
Pmdwn of pressure in the intake pipe with the 
opening of an intake valve corresponding to each 

cylinder, a cylinder filled air amount estimating means for estimating a cylinder filled air 
amount Mti for each cylinder in accordance with the throttle passing air amount detected 
by the throttle passing air amount detecting means and the extra air amount calculated by 
the extra air amount calculating means, and an engine control means for controlling the 
internal combustion engine in accordance with the cylinder filled air amount for each 
cylinder estimated by the cylinder filled air amount estimating means. 
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♦ NOTICES * 

JPO and NCI PI are not: responsible for any 
damages caused by the use of this translation. 

1. This document has been translated by computer. So the translation may not reflect the original 
precisely. 

2. **** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 



CLAIMS 



[Claim(s)] 
[Claim 1] 

A throttle passage air content calculation means to compute the throttle passage air content which 
passes a throttle valve, 

An excessive air content calculation means to compute the excessive air content to the above- 
mentioned gas column equivalent to the amount of descent of the pressure of inhalation of air by the 
inlet valve corresponding to each gas column opening, 

A restoration air content presumption means in a cylinder to presume the restoration air content in a 
cylinder for every gas column based on the throttle passage air content detected by the above- 
mentioned throttle passage air content detection means, and the excessive air content computed by 
the above-mentioned excessive air content calculation means, 

The control imit of the internal combustion engine possessing the engine control means which 
controls ah intemal combustion engine based on the restoration air content in a cylinder for every gas 
column presimied by this restoration air content presumption means in a cylinder. 
[Claim 2] 

The above-mentioned restoration air content presumption means in a cylinder is the control imit of 
the intemal combustion engine according to claim 1 which adopts what totaled the above-mentioned 
throttle passage air content and the excessive air content to each above-mentioned gas column as a 
restoration air content in a cylinder to each gas column. 
[Claim 3] 

The above-mentioned restoration air content presumption means in a cylinder is the control unit of 
the intemal combustion engine according to claim 1 which adopts what continued for two or more 
cycles and equalized what totaled the above-mentioned throttle passage air content and the excessive 
air content to each above-mentioned gas colimm for every gas column as a restoration air content in 
a cylinder to the above-mentioned gas column. 
[Claun 4] 

The pressure sensor for detecting the pressure of inhalation of air is provided. The above-mentioned 
excessive air content calculation means The difference of the maximum of the pressure of inhalation 
of air and the minimum value which were detected by the above-mentioned pressure sensor in the 
period which the inlet valve corresponding to each gas colxmm is opening, and the period of the near. 
The control unit of an intemal combustion engine given in any 1 term of claims 1-3 which compute 
the excessive air content to the above-mentioned gas column using an equation of state based on the 
temperature of inhalation of air. 
[Claim 5] 

The control unit of the intemal combustion engine according to claim 4 which adopts atmospheric 
temperature as the above-mentioned temperature of inhalation of air. 
[Claim 6] 

The above-mentioned excessive air content calculation means is the control unit of an intemal 
combustion engine given in any 1 term of claims 1-3 which compute the excessive air content to the 
above-mentioned gas column based on the amount of rises of the pressure of inhalation of air 
immediately after closing this inlet valve just before the inlet valve corresponding to the amoimt of 
descent and the above-mentioned gas column of the pressure of inhalation of air by the inlet valve 



http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_web_cgi_ejje?u=http%3A%2F%2Fwww4.ip... 9/21/2006 



JP,2005-090437,A [CLAIMS] 



Page 2 of 2 



corresponding to each gas column opening opens. 
[Claim 7] 

The flow rate sensor which detects the throttle-valve passage air flow rate which passes a throttle 
valve is provided. The above-mentioned throttle passage air content calculation means In the period 
during the minimxmi value stage when the maximum stage when the pressure of inhalation of air 
serves as max in the period which the inlet valve corresponding to each gas column is opening, and 
the period of the near, and the pressure of inhalation of air serve as min The control unit of an 
internal combustion engine given in any 1 term of claims 1-6 which compute a throttle passage air 
content by integrating with the throttle-valve passage air flow rate detected by the above-mentioned 
flow rate sensor. 
[Claim 8] 

The flow rate sensor which detects the throttle-valve passage air flow rate which passes a throttle 
valve is provided. The period during the minimum vdue stage when the maximum stage when the 
pressure of inhalation of air serves as max in the period which the inlet valve corresponding to each 
gas column is opening, and the period of the near, and the pressure of inhalation of air serve as min 
is set to deltatdwn. If the period between the valve-opening stage of the above-mentioned inlet valve 
and a clausilium stage is set to deltatioc and the throttle-valve passage air flow rate detected by the 
above-mentioned flow rate sensor during these periods is set to mt The above-mentioned throttle 
passage air content calculation means is the control unit of an intemal combustion engine given in 
any 1 term of claims 1-6 which compute the throttle passage air content Mt based on the following 
type(l). 

Mt=mt- (delta tdwn-f delta tioc)/2 - (1) 
[Claim 9] 

The above-mentioned engine control means is the control unit of an intemal combustion engine 
given in any 1 term of claims 1-8 which control fuel oil consumption and ignition timing for every 
gas colxunn based on the restoration air content in a cylinder for every gas column presiraied by the 
restoration air content presumption means in a cylinder. 
[Claim 10] 

A working angle is changed according to engine operational status, and the above-mentioned inlet 
valve memorizes beforehand the relation between the restoration air content in a cylinder in specific 
engine operational status, and the working angle of the above-mentioned inlet valve. Based on the 
relation by which storage was carried out [ above-mentioned ] with the restoration air content in a 
cylinder computed by the above-mentioned restoration air content calculation means in a cylinder, 
the actual working angle in each gas column is presumed. The control imit of an intemal combustion 
engine given in any 1 term of claims 1-9 which amend an intemal combustion engine's operation 
parameter so that tiie difference of these working angles may be compensated when the presumed 
this actual working angle differs firom a target working angle. 
[Claim 11] 

An air content prediction means to predict the average restoration air content in a cylinder of all gas 
colimms based on throttle opening, and the atmospheric temperature and the atmospheric pressure 
force of the perimeter of an intemal combustion engine at least is provided further. The relative 
deflection between gas columns is computed based on the restoration air content in a cylinder for 
every gas column presumed by the above-mentioned restoration air content presumption means in a 
cylinder when engine operational status was in a steady state. The above-mentioned engine control 
means The control unit of an intemal combustion engine given in any 1 term of claims 1-10 which 
control an intemal combustion engine based on the restoration air content in a cyUnder for every gas 
column computed by having amended the average restoration air content in a cylinder predicted by 
the above-mentioned air content prediction means when engine operational status was in a transient 
based on the above-mentioned deflection. 



[Translation done.] 
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♦ NOTICES * 

JPO smd NCIPI are not responsible for any 
damages caused by the use o£ this translation. 

1 .This document has been translated by computer. So the translation may not reflect the original 
precisely. 

2.**** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 

DETAILED DESCRIPTION 

[Detailed Description of the Invention] 

[Field of the Invention] 

[0001] 

This invention relates to an intemal combustion engine's control unit. 

[Backgroimd of the Invention] 

[0002] 

In order to make into the optimal value the air-fuel ratio of the gaseous mixture which bums in an 
intemal combustion engine's combustion chamber, when an inlet valve closes, it is necessary to 
presxime correctly the amount (for "the restoration air content in a cylinder" to be called hereafter) of 
the air with which the combustion chamber is filled up. Usually, the restoration air content in a 
cylinder is presumed from the map of a large number which made the argument the output value 
from many sensors, such as a flow rate sensor (air flow meter), and these sensors. Here, if the 
restoration air content in a cylinder is presmned using a map, the number of required maps and the 
number of the arguments will increase, and the adaptation man day of map creation time will 
increase very much in connection with this. Then, in recent years, computing the restoration air 
content in a cylinder by reducing the number and argument of a map is examined by using the 
numerical calculation model expressed with the formula based on hydrodynamics etc. 
[0003] 

The equipment which computes the restoration air content in a cyUnder using such a numerical 
calculation model is indicated by the patent reference 1 . The equipment of the patent reference 1 is 
computing the restoration air content in a cylinder using what subtracted the air content stored within 
inhalation of air from the air content which flows into an inlet pipe by the conservation of mass 
being equal to the air content with which it fills up in a cylinder. What subtracted the variation of the 
air of inhalation of air specifically computed from the throttle passage air content detected by the air 
flow meter etc. based on the pressure of inhalation of air detected by the pressTire sensor etc. is 
computed as a restoration air content in a cylinder. 
[0004] 

[Patent reference 1] JP,2002-70633,A 
[Patent reference 2] JP,2001-234798,A 
[Description of the Invention] 
[Problem(s) to be Solved by the Invention] 
[0005] 

By the way, since the inlet valve corresponding to each gas column is opened one by one, pulsation 
occurs to the pressure of inhalation of air (inhalation-of-air pulsation). However, since count will 
become complicated if the restoration air content in a cylinder is computed taking pulsation of the 
pressure of inhalation of air into consideration, he is trying to compute the restoration air content in a 
cylinder by pulsation of the actually produced pressure of inhalation of air being disregarded with 
the equipment of the patent reference 1. Namely, although the pressure of inhalation of air is 
changing with inhalation-of-air pulsation a lot in fact, it is computing the variation of the air of 
inhalation of air by removing change of the pressure of inhalation of air by pulsation by coxmt. 
[0006] 

However, in fact, if pulsation of the pressure of inhalation of air is closely related to the restoration 
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air content in a cylinder to each gas column and the restoration air content in a cylinder is 
computable using such pulsation, it can compute the restoration air content in a cylinder more 
correctly. 
[0007] 

Then, the purpose of this invention presumes the restoration air content in a cylinder for every gas 
column easily using pulsation of the pressure of inhalation of air, and is to offer the control unit of 
the internal combustion engine which can control an intemal combustion engine the optimal. 
[Means for Solving the Problem] 
[0008] 

In order to solve the above-mentioned technical problem, in the 1st invention A throttle passage air 
content calculation means to compute the throttle passage air content which passes a throttle valve. 
An excessive air content calculation means to compute the excessive air content to the above- 
mentioned gas colunm equivalent to the amount of descent of the pressure of inhalation of air by the 
inlet valve corresponding to each gas colunm opening, A restoration air content presimiption means 
in a cylinder to presimie the restoration air content in a cylinder for every gas column based on the 
throttle passage air content detected by the above-mentioned throttle passage air content detection 
means, and the excessive air content computed by the above-mentioned excessive air content 
calculation means. The control unit of the intemal combustion engine possessing the engine control 
means which controls an intemal combustion engine based on the restoration air content in a 
cylinder for every gas column presumed by this restoration air content presumption means in a 
cylinder is offered. 
[0009] 

In the 2nd invention, what totaled the above-mentioned throttle passage air content and the excessive 
air content to each above-mentioned gas column is used for the above-mentioned restoration air 
content presumption means in a cylinder as a restoration air content in a cylinder to each gas column 
in the 1st invention. 
[0010] 

In the 3rd invention, what continued for two or more cycles and equalized what totaled the above- 
mentioned throttle passage air content and the excessive air content to each above-mentioned gas 
column for every gas column is used for the above-mentioned restoration air content presumption 
means in a cylinder in the 1st invention as a restoration air content in a cylinder to the above- 
mentioned gas column. 
[0011] 

In any 3rd one invention, the pressure sensor for detecting the pressure of inhalation of air is 
provided, the 4th invention — the 1- the above-mentioned excessive air content calculation means 
Based on the difference of maximum and the minimxmi value and the temperature of inhalation of air 
of the pressure of inhalation of air detected by the above-mentioned pressure sensor in the period 
which the inlet valve corresponding to each gas column is opening, and the period of the near, the 
excessive air content to the above-mentioned gas column is computed using an equation of state. 
[0012] 

In the 5th invention, atmospheric temperature is adopted as the above-mentioned temperature of 

inhalation of air in the 4th invention. 

[0013] 

the 6th invention — the 1- in any 3rd one invention, the above-mentioned excessive air content 
C£dculation means computes the excessive air content to the above-mentioned gas colimm based on 
the amoimt of rises of the pressure of inhalation of air immediately after closing this inlet valve, just 
before the inlet valve corresponding to the amount of descent and the above-mentioned gas column 
of the pressure of inhalation of air by the inlet valve corresponding to each gas column opening 
opens. 
[0014] 

The flow rate sensor which detects the throttle-valve passage air flow rate which passes a throttle 
valve in any l-6th one invention in the 7th invention is provided. In the period during the minimum 
value stage when the maximimi stage when the pressure of inhalation of air serves as max in the 
period when the inlet valve corresponding to each gas column is opening the above-mentioned 
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throttle passage air content calculation means, and the period of the near, and the pressure of 
inhalation of air serve as min A throttle passage air content is computed by integrating with the 
throttle-valve passage air flow rate detected by the above-mentioned flow rate sensor. 
[0015] 

The flow rate sensor which detects the throttle-valve passage air flow rate which passes a throttle 
valve in any 6th one invention is provided, the 8th invention ~ the 1- The period during the 
minimum value stage when the maximum stage when the pressure of inhalation of air serves as max 
in the period which the inlet valve corresponding to each gas column is opening, and the period of 
the near, and the pressure of inhalation of air serve as min is set to deltatdwn. If the period between 
the valve-opening stage of the above-mentioned inlet valve and a clausilium stage is set to deltatioc 
and the throttle-valve passage air flow rate detected by the above-mentioned flow rate sensor during 
these periods is set to mt The above-mentioned throttle passage air content calculation means is the 
control unit of an internal combustion engine given in any 1 term of claims 1-6 which compute the 
throttle passage air content Mt based on the following type (1). 
Mt=mt- (delta tdwn+delta tioc)/2 - (1) 
[0016] 

In any l-8th one invention, the above-mentioned engine control means controls fuel oil consumption 
and ignition timing by 9th invention for every gas column based on the restoration air content in a 
cylinder for every gas column presumed by the restoration air content presumption means in a 
cylinder. 
[0017] 

In the 10th invention, as for the above-mentioned inlet valve, a working angle is changed in any 1- 
9th one invention according to engine operational status. The relation between the restoration air 
content in a cylinder in specific engine operational status and the working angle of the above- 
mentioned inlet valve is memorized beforehand. Based on the relation by which storage was carried 
out [ above-mentioned ] with the restoration air content in a cylinder computed by the above- 
mentioned restoration air content calculation means in a cylinder, the actual working angle in each 
g£is column is presumed. When the presumed this actual working angle differs from a target working 
angle, an internal combustion engine's operation parameter is amended so that the difference of these 
working angles may be compensated. 
[0018] 

In the 1 1th invention, it sets to any 1-1 0th one invention. At least Throttle opening. An air content 
prediction means to predict the average restoration air content in a cylinder of all gas columns based 
on the atmospheric temperature and the atmospheric pressure force of the perimeter of an intemal 
combustion engine is provided further. The relative deflection between gas columns is computed 
based on the restoration air content in a cylinder for every gas colunm presumed by the above- 
mentioned restoration air content presvmiption means in a cylinder when engine operational status 
was in a steady state. The above-mentioned engine control means When engine operational status is 
in a transient, an intemal combustion engine is controlled based on the restoration air content in a 
cylinder for every gas colimin computed by having amended the average restoration air content in a 
cylinder predicted by the above-mentioned air content prediction means based on the above- 
mentioned deflection. 
[Effect of the Invention] 
[0019] 

According to this invention, a throttle passage air content is computed by the throttle passage air 
content calculation means, and an excessive air content is computed by the excessive air content 
calculation means, the restoration air content in a cylinder is presumed for every gas column based 
on these, and an intemal combustion engine is controlled based on this. The restoration air content in 
a cylinder is computable only from the air content and throttle passage air content equivalent to the 
amount of descent of the pressure of inhalation of air produced by pulsation of the pressure of 
inhalation of air. Therefore, according to this invention, the restoration air content in a cylinder for 
every gas column can be easily presumed using pulsation of the pressure of inhalation of air, and an 
intemal combustion engine can be controlled the optimal. 
[Best Mode of Carrying Out the Invention] 
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[0020] 

Hereafter, the first operation gestalt of this invention is explained with reference to a drawing. The 
engine body 1 roughly shown in drawing 1 shows the injection mold jump-spark-ignition type 
internal combustion engine in a cylinder. However, this invention may be applied to another jump- 
spark-ignition type internal combustion engine or a compression autohesion fire type internal 
combustion engine. 
[0021] 

As shown in drawing 1 , with the first operation gestalt of this invention, the engine body 1 possesses 
a cylinder block 2, the piston 3 which reciprocates within a cylinder block 2, and the cylinder head 4 
fixed on the cylinder block 2. A combustion chamber 5 is formed between a piston 3 and the 
cylinder head 4. An inlet valve 6, a suction port 7, an exhaust valve 8, and the exhaust air port 9 are 
arranged for every gas column at the cylinder head 4. Furthermore, as shown in drawing 1 , an 
ignition plug 10 is arranged in the center section of the internal surface of the cylinder head 4, and a 
fuel injection valve 1 1 is arranged at a cylinder head 4 internal-surface periphery. Moreover, the 
cavity 12 prolonged fi-om the lower part of a fiiel injection valve 1 1 to the lower part of an ignition 
plug 10 is formed in the top face of a piston 3. 
[0022] 

The suction port 7 of each gas column is connected with a surge tank 14 through the inhalation-of-air 
branch pipe 13, and a surge tank 14 is connected with an air cleaner 16 through an inlet pipe 15. In 
an inlet pipe 15, the throttle valve 18 driven with a step motor 17 is arranged. Moreover, the air flow 
meter 19 for detecting the flow rate of the air (inhalation-of-air gas) which passes an inlet pipe 15 to 
the inlet pipe 15 of the throttle- valve 18 upstream is arranged. On the other hand, the exhaust air port 
9 of each gas column is connected with an exhaust pipe 20, and this exhaust pipe 20 is connected 
with an exhaust emission control device 21 . 
[0023] 

An electronic control unit (ECU) 31 consists of a digital computer, and RAM (random access 
memory)33, ROM (read-only memory)34, CPU (microprocessor)35, the input port 36, and the 
output port 37 which were mutually connected through the bidirectional bus 32 are provided. The 
temperature sensor 41 of inhalation of air for detecting the temperature of the air the pressure sensor 
40 of inhalation of air for detecting the pressure of the air within inhalation of air (inhalation-of-air 
gas) and within inhalation of air in a surge tank 14 is formed. The pressure sensor 40 of these 
inhalation of air and the temperature sensor 41 of inhalation of air generate the output voltage which 
is proportional to the pressure of inhalation of air, and the temperature of inhalation of air, 
respectively, and are inputted into input port 36 through A-D converter 38 to which this output 
voltage corresponds. 
[0024] 

Moreover, the throttle opening sensor 42 for detecting the opening of a throttle valve 1 8, The 
atmospheric temperature sensor 43 for detecting the atmospheric temperature around an intemal 
combustion engine, or the temperature (intake-air temperature) of the air inhaled by the inlet pipe 15, 
The atmospheric pressure sensor 44 for detecting the atmospheric pressure force aroimd an intemal 
combustion engine or the pressure (intake pressure) of the air inhaled by the inlet pipe 15 is formed, 
and the output voltage of these sensors is inputted into input port 36 through corresponding A-D 
converter 38. Moreover, the load sensor 46 which generates flie output voltage proportional to the 
amount of treading in of an accelerator pedal 45 is connected to an accelerator pedal 45, and the 
output voltage of the load sensor 46 is inputted into input port 36 through corresponding A-D 
converter 38. The crank angle sensor 47 generates an output pulse, whenever a crankshaft rotates 30 
degrees, and this output pulse is inputted into input port 36. In CPU35, an engine rotational 
frequency is calculated from the output pulse of this crank angle sensor 47. On the other hand, an 
output port 37 is connected to an ignition plug 10, a fiiel injection valve 11, and a step motor 17 
through the corresponding drive circuit 39. 
[0025] 

by the way, in order to make into a target air-ftiel ratio the air-fiiel ratio of the gaseous mixture which 
bums in an intemal combustion engine's combustion chamber 5 in an intemal combustion engine's 
control unit The amount of the air (inhalation-of-air gas) with which it fills up in the combustion 
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chamber 5 when an inlet valve closes The amoxrnt ("fuel oil consumption" is called hereafter) of the 
fuel injected from a fiiel injection valve to an internal combustion engine's combustion chamber 5 (or 
inhalation-of-air path) so that it may presume ("the restoration air content Mc in a cylinder" is called 
hereafter) and the air-fiiel ratio of gaseous mixture may turn into a target air-fiiel ratio based on the 
presumed restoration air content Mc in a cylinder is defined. Therefore, in order to make correctly 
into a target air-fiiel ratio the air-fiiel ratio of the gaseous mixture which bums in an internal 
combustion engine's combustion chamber S, it is necessary to presume the restoration air content Mc 
in a cylinder correctly. 
[0026] 

Usually, the restoration air content Mc in a cylinder is presumed from the map of a large number 
which made the argument the output value from many sensors, such as a flow rate sensor (air flow 
meter), and these sensors. However, when presuming the restoration air content Mc in a cylinder 
using a map in this way, in order to make more exact the value of the restoration air content Mc in a 
cylinder presumed, the number of required maps and tiie number of the arguments increase. Thus, if 
the number of maps increases, ROM of ECU for saving a map will have to be made into what has 
large storage capacity, and the manufacturing cost of an internal combustion engine's control unit 
will become high. Furthermore, although an adaptation activity must be done for every format of the 
internal combustion engine by which a map is used for creating each map, since the point of 
measurement in this adaptation activity increases according to the number of maps, and the number 
of those arguments, if the number of maps and the number of those arguments increase, the man day 
of an adaptation activity will also increase. 
[0027] 

Then, the control unit of the intemal combustion engine which computes the restoration air content 
Mc in a cylinder by numerical calculation using various models, without using a map is examined. In 
such a control unit, though the man day at the time of trying to reduce the number of required maps 
as much as possible, and this doing an adaptation activity by using numerical calculation abundantly 
is reduced sharply, the restoration air content Mc in a cylinder is correctly computable. 
[0028] 

Some which compute the restoration air content Mc in a cylinder are in one of such the models from 
the flow rate ("the throttle-valve passage air flow rate mt" is called hereafter) of the air which passes 
a throttle valve 1 8 to per unit time amount, and the pressure ("the pressure Pm of inhalation of air" is 
called hereafter) of the air which exists in the part ("a part for an inhalation-of-air tube part" is called 
hereafter) of the inlet-pipe 15 grade from a throttle valve 18 to an inlet valve 6 (for example, patent 
reference 1). The flow rate of the air inhaled in a cylinder in such a model ("the intake air flow mc in 
a cylinder" is called hereafter) 0 That is, the flow rate of the inhalation-of-air gas which flows out of 
a part for an inhalation-of-air tube part is the throttle- valve passage air flow rate mt (). That is, the 
conservation of mass of being equal to what subtracted the amount (namely, amount of the 
inhalation-of-air gas stored in the inhalation-of-air tube part part) of the inhalation-of-air gas 
equivalent to the amoxrnt of rises of the pressure Pm of inhalation of air in unit time amount from the 
flow rate of the air which flows into a part for an inhalation-of-air tube part is used. 
[0029] 

Usually, the inhalation-of-air pulsation by an inlet valve being opened one by one in a part for an 
inhalation-of-air tube part occurs, and, tiierefore, the pressure of inhalation of air is changed sharply. 
Thus, if the model using the conservation of mass mentioned above using the pressure of inhalation 
of air changed sharply is created, a model type will become complicated and a count load will 
become large. So, in the former, in order to eliminate the effect of fluctuation of the pressure of 
inhalation of air by inhalation-of-air pulsation, the detection value and detection value of for 
example, the pressiu*e sensor of inhalation of air aimeal as variation (dPm/dt) of the pressure of 
inhalation of air per unit time amount, and deflection with a value is used. 
[0030] 

However, if the restoration air content in a cyUnder is computed by fluctuation of the pressure of 
inhalation of air by inhalation-of-air pulsation having influenced the restoration air content in a 
cylinder greatly, therefore disregarding this effect, it caimot compute the exact restoration air content 
in a cylinder. Conversely, fluctuation of the pressure of inhalation of air by inhalation-of-air 
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pulsation if it says can compute correctly the restoration air content in a cylinder to each gas column, 
if it uses that it is closely related to the restoration air content in a cylinder. So, in this invention, the 
restoration air content in a cylinder is computed using this. 
[0031] 

Hereafter, with reference to drawing 2 and drawing 3 , the calculation approach of the restoration air 
content in a cyUnder is explained. In addition, drawing 2 shows the fundamental concept of the 
model (a "inlet-pipe model" is called hereafter) Ml in a part for an inhalation-of-air tube part. 
Drawing 3 (a) shows change of the flow rate to a crank angle. The continuous line mt in drawing 3 
shows the throttle- valve passage air flow rate, and the continuous line mci shows the intake air flow 
in a cylinder to all gas columns. Moreover, drawdng 3 (b) shows change of the pressure of inhalation 
of air to a crank angle. 
[0032] 

First, the inlet-pipe model Ml shown in drawing 2 is considered. If the conservation of mass is 
applied about a part for an inhalation-of-air tube part, the relation of the following formula (2) will 
be realized to the pressure Pm of inhalation of air, the air flow rate (namely, throttle-valve passage 
air flow rate mt) which flows into a part for an inhalation-of-air tube part, and the flow rate (namely, 
intake air flow mci in a cylinder to a No. i gas colimm) of the inhalation-of-air gas which flows out 
of a part for an inhalation-of-air tube part. 
[Equation 1] 



dPm Ra Tm . ^ ^ 

{mt—y^mci) •••(2) 



dt Vm 



It is the value to which the temperature ("the temperature of inhalation of air" is called hereaifter) of 
the air to which Tm exists in an inhalation-of-air tube part part, and Vm did the gas constant to Ihe 
volume for an inhalation-of-air tube part here, and Ra did the division with the average molecular 
weight of air. Therefore, variation deltaPm of the pressure of inhalation of air for [ delta ] t seconds 
can be expressed like the following formula (3) by integrating with a formula (2) from time of day t. 
[Equation 2] 



APm = I (mt -2\mci) (3) 

Vm ^' ^ 



[0033] 

If the inflow air content (mt) for an inhalation-of-air tube part is larger than an outflow air flow rate 
(mci), the pressure of inhalation of air will rise, if small, the pressure of inhalation of air will 
descend, if equal, the pressure of inhalation of air is fixed, and a formula (3) shows that variation 
deltaPm of the pressure of inhalation of air for [ delta ] t seconds is equivalent to a changed part of 
the air content in an inhalation-of-air tube part part. In addition, when it is in a steady state which 
engine operational status mentions later, the amount of inhalation-of-air tube part is buffered to 
being intermittent depending on closing motion of an inlet valve 6, and the change is [ the inflow air 
flow rate (mt) from a part for an inhalation-of-air tube part ] quiet [ the outflow air flow rate (mci) 
from a part for an inhalation-of-air tube part ]. For this reason, the size relation between an outflow 
air flow rate (mci) and an inflow air flow rate (mt) repeats reversal (refer to drawing 3 (a)). This 
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means that the value in the parenthesis of the right-hand side of the above-mentioned fomiula (2) 
repeats positive/negative reversal a fixed period, i.e., the pressure of inhalation of air repeats a rise 
and descent a fixed period, and expresses pulsation of the pressure of inhalation of air. 
[0034] 

Here, it is assumed that the valve-opening period of the inlet valve 6 of each gas column does not 
overlap as shown in drawing 3 (a). That in this case, the pressure of inhalation of air takes Maximum 
Pmmax about the inhalation to tiie i-th cylinder When the time amount differential value of the 
pressure of inhalation of air is zero (dPm/dt=0), it is a time (mt=mci) of the magnitude of the 
throttle-valve passage air flow rate mt and the cylinder [ i-th ] intake air flow mci in a cylinder 
balancing. It is at the time whose throttle- valve passage air flow rate mt was larger while the amount 
mci of inspired gas in a cylinder is increasing, namely, until the above-mentioned magnitude 
balances (let time of day be the maximimi time of day tmax for this time). On the other hand, it is a 
time of the time amount differential value of the pressure of inhalation of air being zero that the 
pressure of inhalation of air takes the minimima value Pmmin about the inhalation to the i-th 
cylinder, and it is at the time whose amount mci of inspired gas in a cylinder was larger while the 
amount mci of inspired gas in a cylinder is decreasing, namely, imtil ttie above-mentioned magnitude 
balances (let time of day at this time be the minimum value time of day tmin). 
[0035] 

Therefore, amount of descent ("amount of pressure drawdowns of inhalation of air" is called 
hereafter) deltaPmdwn (namely, difference of the maximxmi Pmmax of the pressure of inhalation of 
air and the minimum value Pmmin) of the pressure of inhalation of air produced by inhalation of 
cylinder [ i-th ] inhalation-of-air gas can be expressed like the following type (4). In addition, the 
integral term of a formula (4) is equivalent to tiie area A of drawing 3 (a), and it turns out that 
deltaPmdwn is proportional to area A. Therefore, the capacity corresponding to area A can be called 
the excessive capacity to the i-th cylinder equivalent to the amount of descent of the pressure of 
inhalation of air by the inlet valve corresponding to the i-th cylinder opening. 
[Equation 3] 



APmdwn = 



Ra-Tm r""^ 



Vm 



(mt - y mci)dt 



(4) 



[0036] 

From an assumption that the valve-opening period of the inlet valve 6 of each gas column does not 
overlap, the above-mentioned formula (4) can deform like the following formula (5). 
[Equation 4] 



Jtmax / I j Jtmax 



Here, the integral term of the throttle- valve passage air flow rate mt in a formula (5) is equivalent to 
the area B of drawing 3 (a), and Mci serves as a value adding the area A and area B of drawing 3 . 
Therefore, Mci is equivalent to the capacity with which it filled up in the cylinder [ i-th ] combustion 
chamber 5 during the valve-opening period of the inlet valve 6 corresponding to the i-th cylinder, 
i.e., the restoration air content in a cylinder. However, if it says strictly, since the actual restoration 
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air content in a cylinder is equivalent to the amount which applied area C to the area A and area B of 
drawing 3 , Above Mci serves as approximate value which disregarded the capacity equivalent to 
area C as very small 
[0037] 

Therefore, when engine operational status is in a steady state and the valve-opening period of the 
inlet valve 6 of each gas column does not overlap, the cyhnder [ i-th ] restoration air content Mci in a 
cylinder can be presumed from the above-mentioned formula (5) by detecting or computing the 
throttle- valve passage air flow rate mt, the temperature Tm of inhalation of air, and amount of 
pressure drawdowns of inhalation of air deltaPmdwn. 
[0038] 

In addition, in mounting a formula (5), a formula (5) may be transformed like the following type (6). 
[Equation 5] 



Mci = tsJPmdwn + • Ltdwn • - (6) 

/ K Vm J 



In a formula (6), deltatdwn is the time amount of a before [ from the maximum time of day tmax / 
the minimum value time of day train ], and expresses the decay time of the pressure of inhalation of 
air. Moreover, the throttle-valve passage air flow rate mt in a formula (6) is the value which 
averaged the detection value of an air flow meter 19 during the period from the maximum time of 
day tmax to the minimum value time of day tmin, or the valve-opening period of an inlet valve 6. Or 
since fluctuation of the actual throttle-valve passage air flow rate in the period from the maximum 
time of day tmax to the minimum value time of day tmin or the valve-opening period of an inlet 
valve 6 is small, you may be the detection value of the air flow meter 19 in the specific time of day 
in these periods. Similarly, the temperature Tm of inhalation of air in a formula (6) is also the value 
which averaged the detection value of the temperature sensor 41 of inhalation of air in the above- 
mentioned period, or the detection value of the temperature sensor 41 of inhalation of air in the 
specific time of day in the above-mentioned period. 
[0039] 

Moreover, although the temperature sensor 41 of inhalation of air was attached in the surge tank 14 
and the temperature of the inhalation-of-air gas in an inhalation-of-air tube part part is detected with 
the above-mentioned operation gestalt, a temperature sensor may be attached in the inhalation-of-air 
upstream of a throttle valve 18, or a temperature sensor may be formed in one with an air flow meter 
19, and the temperature detected by this temperature sensor may be used as temperature of inhalation 
of air. It is because it is carried out when especially this can approximate [ temperature / of the air of 
the inhalation-of-air upstream of a throttle valve 18 ] the temperature of inhalation of air with an 
equal mostly when engine operational status is in a steady state, and presumption of the cylinder [ i- 
th ] restoration air content Mci in a cylinder has engine operational status in a steady state with this 
operation gestalt. 
[0040] 

With reference to drawing 4 , the operating procedure of presumption of the cylinder [ using the 
formula (5) of the inlet-pipe model Ml mentioned above / i-th ] restoration air content in a cylinder 
is explained. In addition, as for this actuation, it is desirable to perform, when it performs for every 
predetermined time interval and every gas column, especially engine operational status is in a steady 
state and the valve-opening period of the inlet valve 6 of each gas column does not overlap. 
[0041] 

First, at step 101, 1 is added to the time amount counter n. After a cylinder [ i-th ] inlet valve closes 
the valve last time, the count of activation of this actuation is expressed, therefore the time amoimt 
coimter n expresses the elapsed time from an inhalation-of-air valve-closing valve. Hereafter, the 
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value of a time amount counter is explained as time of day. Subsequently, at step 102, current crank 
angle CA is acquired from the crank angle sensor 47. At step 103, the value set to the closing motion 
valve flag Vlv (n) from crank angle CA acquired in step 102 is computed. In addition, the closing 
motion valve flag Vlv (n) expresses the closing motion valvate form voice of the cylinder [ i-th ] 
inlet valve 6 in time of day n, while the cylinder [ i-th ] inlet valve 6 is opening at tiiis time, the value 
of the closing motion valve flag Vlv (n) is set to 1, and while closing the valve, it is set to 0. 
[0042] 

Subsequently, at step 104, it is judged whether it is in the condition that the value of the closing 
motion valve flag Vlv (n-1) in time of day n-1 is set to 1, and the value of the closing motion valve 
flag Vlv (n-1) in time of day n is set to 0. That is, at step 104, it is judged whether whether the inlet 
valve's 6 opening at the time of the last actuation, and the inlet valve's 6 closing the valve at the time 
of this actuation and the time of this actuation that is, are [ inlet valves 6 ] the times of clausiliimi 
being carried out. In step 104, the time of this actuation progresses to step 105, when it judges that an 
inlet valve 6 is not a time of clausilium being carried out. 
[0043] 

At step 105, it is judged whether the inlet valve (whether the value of the closing motion valve flag 
Vlv in time of day n (n) is 0 and the i-th cylinder) 6 is opened. When judged with the cylinder [ i-th ] 
inlet valve 6 not being opened (Vlv(n) =0), step 106 - step 1 12 are not performed, but this actuation 
is made to end them. 
[0044]. 

On the other hand, when judged with the cylinder [ i-th ] inlet valve 6 being opened in step 105, it 
progresses to step 106. At step 106, 1 is added to the valve-opening counter m. In addition, after an 
inlet valve 6 is opened, the count of activation of this actuation is expressed, therefore the valve- 
opening counter m expresses the elapsed time from inlet- valve 6 valve opening. At step 107, the 
pressure Pm of inhalation of air, the throttle-valve passage air flow rate mt, and the temperature Tm 
of inhalation of air are acquired from the pressure sensor 40 of inhalation of air, an air flow meter 19, 
and the temperature sensor 41 of inhalation of air, respectively. 
[0045] 

At steps 108-1 1 1, the maximum Pmmax of the pressure of inhalation of air in an inlet-valve 6 valve- 
opening period and the minimum value Pmmin, the maximum stage tmax, and the minimimi value 
stage tmin are updated. 

Step 109 is performed only when judged with the pressure Pm of inhalation of air which it was 
judged whether whether the pressure Pm of inhalation of air acquired at step 107 being larger than 
the maximum Pmmax of the pressure of inhalation of air by which current storage is carried out, and 
the acquired pressure Pm of inhalation of air are maxes after an inlet valve 6 opens, and it acquired 
being max at step 108 (Pm>Pmmax). At step 109, the pressure Pm of inhalation of air acquired at 
step 107 is memorized as maximum Pmmax of the pressure of inhalation of air, and the current time 
of day n is memorized as maximum time of day tmax. 
[0046] 

Subsequently, step 1 1 i is performed only when judged with the pressure Pm of inhalation of air 
which it was judged whether whether the pressure Pm of inhalation of air acquired at step 107 being 
smaller than the minimum value Pmmin of the pressure of inhalation of air by which current storage 
is carried out, and the acquired pressure Pm of inhalation of air are mins after an inlet valve 6 opens, 
and it acquired at step 110 being min (Pm<Pmmin). At step 1 1 1 , the pressure Pm of inhalation of air 
acquired at step 107 is memorized as the minimum value Pmmin of the pressure of inhalation of air, 
and the current time of day n is memorized £ts minimiun value time of day tmin. 
[0047] 

At step 112, this throttle- valve passage air flow rate acquired at step 107 is added to addition value 
sigmamt of the throttle-valve passage air flow rate after the cylinder [ i-th ] inlet valve 6 opens. 
Furthermore, this temperature Tm of inhalation of air acquired at step 107 is added to addition value 
sigmaTm of the temperature of inhalation of air after the cylinder [ i-th ] inlet valve 6 opens. 
[0048] 

On the other hand, when judged with it being a time of clausilium of the inlet valve 6 being carried 
out for the time of this actuation in step 104, it progresses to step 113. Let the value which did the 
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division of the addition value sigmamt of the throttle- valve passage air flow rate in an inlet- valve 6 
valve-opening period with the value m of a valve-opening counter be the average throttle-valve 
passage air flow rate mtave at step 113. This average throttle- valve passage air flow rate mtave 
expresses the average of the throttle-valve passage air flow rate in the valve-opening period of an 
inlet valve 6. Moreover, the value which did the division of the addition value sigmaTm of the 
temperature of inhalation of air in an inlet-valve 6 valve-opening period with the value m of a valve- 
opening counter is made into the temperature Tmave of average inhalation of air. This temperature 
Tmave of average inhalation of air expresses the average of the temperature of inhalation of air in 
the valve-opening period of an inlet valve 6. 
[0049] 

Subsequently, at step 1 14, the value which subtracted the minimum value Pmmin of the pressure of 
inhalation of air updated in step 111 from the maximum Pmmax of the pressure of inhalation of air 
updated in step 109 is set to amount of pressure drawdowns of inhalation of air deltaPmdwn 
(deltaPmdwn=Pmmax-Pmmin). At step 1 15, the value which subtracted the maximimi time of day 
tmax updated in step 109 from the minimum value time of day tmin updated in step 1 1 1 is set to 
deltatdwn (deltatdwn=tmin-tmax), 
[0050] 

At step 116, the restoration air content Mci in a cylinder to the cylinder [ i-th ] combustion chamber 
5 is computed by substituting for a formula (6) mtave and Tmave which were computed at steps 113- 
115, deltaPmdwn, and deltatdwn. Subsequently, at step 117, the value of Coxmters n and m is reset 
by zero, the value of Pnmiin is reset by infinity and the value of addition value sigmamt and 
sigmaTm is reset for the value of Pmmax by zero at zero, respectively. 
[0051] 

By the way, in the control unit of the internal combustion engine of this operation gestalt, the fuel oil 
consumption from the fiiel injection valve 1 1 injected into the i-th cylinder so that the air- fuel ratio 
of the gaseous mixture in the i-th cylinder may turn into a target air- fuel ratio is defined based on the 
cylinder [ which was presumed as mentioned above / i-th ] restoration air content Mci in a cylinder. 
A target air-fiiel ratio is determined by ECU3 1 based on engine operational status (for example, an 
engine rotational frequency and an engine load) etc. thereby, in agreement [ in the air- fuel ratio of 
gaseous mixture / almost correctly ] with a target air-fuel ratio about all gas columns, even if 
dispersion in the restoration air content in a cylinder arises between gas columns ~ it can make — 
exhaust air — aggravation of description can be controlled. 
[0052] 

However, if dispersion in the restoration air content in a cylinder arises between gas colimms when 
fuel oil consumption is determined in this way, the fuel oil consumption injected from a fuel 
injection valve 1 1 between gas columns will differ. For this reason, the combustion energy 
("combustion energy" is only called hereafter) which contributes to generating by combustion of a 
fuel and depressing a piston 3 will also differ between gas columns, and, therefore, torque fluctuation 
will produce it. Therefore, in order to control generating of torque fluctuation, it is necessary to make 
it combustion energy become equal between gas columns in addition to determining fuel oil 
consumption that the air- fuel ratio of the gaseous mixture in each gas column will tmn into a target 
air- fuel ratio as mentioned above. 
[0053] 

So, with this operation gestalt, combustion energy is arranged with homogeneity between gas 
columns by adjusting the ignition timing by the ignition plug 10 for every gas colimm. This situation 
is explained taking the case of the 1st cylinder and the 2nd cylinder with reference to drawing 5 . 
Drawing 5 shows the relation of the ignition timing and combustion energy in each gas column. 
TDC shows the compression top dead center in each gas column among drawing. 
[0054] 

The case where dispersion in the restoration air content in a cylinder is between the 1st cylinder and 
the 2nd cylinder, there are few cylinder [ 1st ] restoration air contents in a cylinder than the cylinder 
[ 2nd ] restoration air content in a cylinder, the relation of the ignition timing and combustion energy 
in the 1st cylinder becomes like continuous-line #1 as shown in drawing 5 , and the relation of the 
ignition timing and combustion energy in the 2nd cylinder has become like continuous-line #2 is 
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considered. 
[0055] 

In this case, it is the target ignition timing (CAl in drawing.) most computed by the ignition timing 
31 by the side of a tooth lead angle, i.e., ECU, in the range which is not generated, such as knocking, 
about the 1st cylinder. Hereafter, ignition according to an ignition plug 10 for calling "the 1st 
ignition timing" is performed. If ignition is perfomied about the 2nd cylinder at the same ignition 
timing as the 1st ignition timing CA 1 at this time, the combustion energy to generate will become 
larger than a cylinder [ 1st ] thing. Therefore, it is the ignition timing when the almost same 
combustion energy as the combustion energy generated in the 2nd cylinder when ignition is 
performed in the 1st cylinder at the 1st ignition timing CA 1 occurs, and is the ignition timing by the 
side of a lag (CA2 in drawing.) from the 1st ignition timing CA 1. Hereafter, ignition is performed 
for calling it "the 2nd ignition timing." the combustion energy generated in each gas column by 
carrying out like this, maintaining the air-fuel ratio of gaseous mixture at homogeneity mostly 
between gas columns — between gas columns — almost — homogeneity — it can maintain — therefore 
— exhaust air — torque fluctuation can be controlled, controlling aggravation of description. 
[0056] 

in addition, although only the 1st cylinder and the 2nd cylinder were made into the example and the 
above-mentioned explanation explained it, it is carried out about all gas columns (this operation 
gestalt — like — the case of a 4-cylinder — all 4-cylinders). Therefore, the ignition timing of other gas 
columns is set as the combustion energy which ignition timing of a gas column with few restoration 
air contents in a cylinder is made target ignition timing among all gas columns, and is generated in 
the gas column becoming equal to the combustion energy generated in other gas columns. 
[0057] 

With reference to drawing 6 , the operating procedure which sets the fuel oil consumption and 
ignition timing in each gas colxmm is explained. In addition, this actuation is performed for every gas 
column and every cycle. First, in step 121, the target air-fuel ratio AFt and the target ignition timing 
CAinjt which were computed by ECU31 based on engine operational status, such as an engine 
rotational frequency and an engine load, are acquired. Subsequently, at step 122, the cylinder [ which 
was computed by the actuation shown in drawing 4 / i-th ] restoration air content Mci in a cylinder is 
acquired. 
[0058] 

At step 123, after acquiring the cylinder [ i-th ] restoration air content Mci in a cylinder in a front 
cycle, the smallest restoration air content Mcmin in a cylinder ("the restoration air content in the 
minimum cylinder" is called hereafter) is computed about all gas columns. For example, when there 
are few cylinder [ 1st ] restoration air contents Mci in a cylinder than the restoration air content in a 
cylinder to all other gas columns, the restoration air content Mcmin in the minimum cylinder turns 
into the cylinder [ 1st ] restoration air content Mci in a cylinder. Subsequently, at step 124, the value 
which did the division of the cylinder [ which was acquired in step 122 / i-th ] restoration air content 
Mci in a cylinder with the target air- fuel ratio AFt acquired in step 121 is made into the fuel oil 
consumption TAUi in the i-th cylinder (TAUi=Mci/AFt), and the fuel of this fiiel oil consumption 
TAUi is injected from the cylinder [ i-th ] fuel injection valve 11 at the time of fuel injection. 
[0059] 

At step 125, amount of lags deltaCAinji of the ignition timing in the i-th cylinder is computed based 
on the difference which subtracted the restoration air content Mcmin in the minimum cylinder 
acquired at step 123 from the cylinder [ which was acquired at step 122 / i-th ] restoration air content 
Mci in a cylinder, in addition, amount of lags deltaCAinji and the above — the relation with 
difference is beforehand computed by count experimentally, it is saved at ROM34 of ECU3 1, and 
this map is used in calculation of amount of lags deltaCAinji in step 125. the case where the cylinder 
[ i-th ] restoration air content Mci in a cylinder is the restoration air content Mcmin in the minimum 
cylinder — the above ~ difference is zero and also let amount of lags deltaCAinji of the ignition 
timing in the i-th cylinder be zero. Subsequently, let the value to which amount of lags deltaCAinji 
of the ignition timing in the i-th cylinder computed at step 125 at the target ignition timing CAinjt 
was added be the cylinder [ i-th ] ignition timing CAinji at step 126 (CAinji=CAinjt+delta CAinji). 
Thus, the cylinder [ i-th ] ignition plug 10 is lit at the computed cylinder [ i-th ] ignition timing 
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CAinji. 
[0060] 

Next, the control unit of the second operation gestalt of this invention is explained. Although the 
control unit of the second operation gestalt is the same as the control unit of the first operation gestalt 
fundamentally, the operating procedure of presumption of the cylinder [ i-th ] restoration air content 
in a cylinder differs. 
[0061] 

By the way, since the noise is contained in the output of the pressure sensor 40 of inhalation of air, 
an error may arise in the value of amount of pressure drawdowns of inhalation of air deltaPdwn 
computed based on the output of the pressure sensor 40 of inhalation of air. In connection with this, 
an error may arise also in the value of the restoration air content Mci in a cylinder computed using 
the value of this deltaPdwn. If fuel oil consumption etc. is defined based on the restoration air 
content Mci in a c^hnder including such an error, the actual air-fuel ratio of gaseous mixture will 
stop being in agreement with a target air-fiiel ratio. 
[0062] 

So, with this operation gestalt, the error of the restoration air content Mci in a cylinder which was 
mentioned above is corrected by taking the value ("the restoration air content Mciave in an average 
cylinder" being called hereafter) which continued for two or more cycles and equalized the 
restoration air content in a cylinder computed by actuation of drawing 4 about each gas column. And 
fuel oil consumption etc. is defined based on this restoration air content Mciave in an average 
cylinder. Even if an error arises in the value of deltaPmdwn by this by containing a noise in the 
output of the pressure sensor 40 of inhalation of air etc., effect of with error to the restoration air 
content in a cylinder presimied can be made small, and, therefore, the actual air-fuel ratio of geiseous 
mixture can be made mostly in agreement with a target air-fuel ratio. 
[0063] 

With reference to drawing 7 , the operating procedure which equalizes and presumes the cylinder [ i- 
th ] restoration air content in a cylinder between cycles is explained. In addition, since it is the same 
as that of step 101 of drawing 4 - step 115, and step 117 respectively, step 141 - step 155, and step 
159 omit explanation. 
[0064] 

At step 156, 1 is added to the value eye of a cycle counter. A cycle counter is a coxmter showing the 
number of cycles from an engine start up. Subsequently, at step 157, the restoration air content Mci 
(eye) in a cylinder in this cycle eye is computed by the above-mentioned formula (5) like step 1 1 6 of 
drawing 4 . 
[0065] 

Subsequently, at step 158, as shown in the following formula (7), the value to which only the 
predetermined number Nave did the division of the sum total of the restoration air content Mci in a 
cylinder from the front cycle (cyc-Nave) to this cycle eye with the predetermined number Nave from 
this cycle eye is computed as a restoration air content Mciave in an average cylinder. 
[Equation 6] 



eye 



Mciave = ^ Mci(k) I Nave • ■ • (7) 

k^cyc-Nave 



Here, a predetermined number Nave is the value defined beforehand. With this operation gestalt, this 
restoration air content Mciave in an average cylinder is used in step 122 of drawing 6 , and is used 
for computing the fuel oil consumption of the i-th cylinder, and ignition timing. 
[0066] 
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In addition, with the second operation gestalt, actuation shown in drawing 7 instead of the actuation 
shown in drawing 4 is performed, and actuation in the first operation gestalt and same actuation are 
performed except it. Moreover, although the restoration air content Mciave in an average cylinder is 
made into the average of the restoration air content Mci in a cylinder of a predetermined number 
Nave with the above-mentioned operation gestalt, it is good also as other values, such as a weighted 
average efficiency. 
[0067] 

Next, the control unit of the internal combustion engine of the third operation gestalt of this 
invention is explained. Although the control device of the third operation gestalt is the same as the 
control device of the first operation gestalt fimdamentally, the temperature sensor 41 of inhalation of 
air is not attached in the inlet pipe containing a surge tank 14 of the inhalation-of-air upstream of a 
part for an inhalation-of-air tube part, or a throttle valve 18. Hereafter, with reference to drawing 8 , 
the presumed approach of the restoration air content in a cylinder in the third operation gestalt is 
explained. In addition, drawing 8 is the same drawing as drawing 3 . 
[0068] 

By the way, if a throttle-valve passage air flow rate when the pressure of inhalation of air is 
descending is set to mtdwn, the above-mentioned formula (6) can be expressed like the following 
type (8). 
[Equation 7] 



Mci = APmdwn/ \ — p j-^- mtdwn * Atdwn •••(8) 



Here, as the gas column which an inlet valve 6 opens cylinder [ i-th ] ago is made into the h-th 
cylinder and it was shown in drawing 8 , after the pressure of inhalation of air takes the minimimi 
value about the inhalation to the h-th cylinder, a period until the pressure of inhalation of air takes 
maximum about the inhalation to the i-th cylinder is defined as deltatup, and the amount of rises of 
the pressure of inhalation of air in this period is defined as deltaPmup. 
[0069] 

Thus, when deltaPmup is defined, the restoration air content in a cylinder (the h-th cylinder in the 
period when the pressure of inhalation of air is rising, and the i-th cylinder) can be approximated if 
equal to about 0. For this reason, the above-mentioned formula (3) can deform like titie following 
type (9), and if a throttle- valve passage air flow rate when the pressure of inhalation of air is rising is 
set to mtup, it can deform like the following formula (10). 
[Equation 8] 



tsPmup = I mt dt (9) 



APm up = miup • Atup •'■(10) 

Vm 



If this formula (10) is transformed and it substitutes for a formxila (8), the following formula (11) can 
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be obtained. 
[Equation 9] 



Mci = mtup • blup ^^^^ + mtdwn • l^dwn — (1 1) 

dkPmup 



[0070] 

Namely, amount of descent deltaPmdwn of the pressure of inhalation of air by opening the cylinder 
[ i-th ] inlet valve 6 according to the formula (11), The throttle- valve passage air flow rate mtdwn 
when decay time deltatdwn of the pressure of inhalation of air and the pressure of inhalation of air 
are descending The cylinder [ i-th ] restoration air content Mci in a cylinder is computable from the 
throttie-valye passage air flow rate mtup when amount of rises deltaPmup of the pressure of 
inhalation of air before opening the cylinder [ i-th ] inlet valve 6, rise-time deltatup of the pressure of 
inhalation of air, and the pressure of inhalation of air are rising. 
[0071] . 

Therefore, according to the third operation gestalt, when the valve-opening stage of the inlet valve 6 
of each gas column does not overlap, it sets. By detecting and computing deltaPmup and deltatup by 
the detection / calculation approach of deltaPmdwn in the first operation gestalt, and deltatdwn, and 
the same approach The restoration air content in a cylinder to each gais column can be computed 
without using what kind of temperature sensor, and, therefore, reduction of a manufacturing cost can 
be aimed at. 
[0072] 

In addition, although [ the above-mentioned operation gestalt ] the cylinder [ i-th ] restoration air 
content Mci in a cylinder is computed based on amount of descent deltaPmdwn of the pressure of 
inhalation of air by opening the cylinder [ i-th ] inlet valve 6, and amount of rises deltaPmup of the 
pressure of inhalation of air before opening the cylinder [ i-th ] inlet valve 6 Instead of amount of 
rises deltaPmup of the pressure of inhalation of air before opening the cylinder [ i-th ] inlet valve 6, 
you may compute based on amount of rises deltaPmup of the inhalation-of-air guidance pressure 
after opening the cylinder [ i-th ] inlet valve 6. 
[0073] 

Next, the control unit of the internal combustion engine of the fourth operation gestalt of this 
invention is explained. The control xmit of ttie first operation gestalt is used about the case where the 
valve-opening stage of the inlet valve 6 of each gas column does not overlap fimdamentally. 
However, if the control unit of the first operation gestalt is used about the case where the valve- 
opening stage of the inlet valve 6 of each gas column overlaps, the error of the restoration air content 
Mci in a cylinder to each gas colunm computed will become large. 
[0074] 

That is, as explained using drawing 3 (a), in the first operation gestalt, the restoration air content in a 
cylinder serves as approximate value which disregarded the edacity equivalent to area C as very 
small. However, when the valve-opening stage of the inlet valve 6 between gas colxmms overlaps, as 
shown in drawing 9 , the throttle- valve passage air flow rate mt is large, and it is so large that the 
capacity which is therefore equivalent to area C cannot be disregarded. 
[0075] 

So, with the fourth operation gestalt, it is supposed that any capacity other than the capacity which is 
equivalent to area A among the restoration air contents Mci in a cylinder to each gas column is not 
calculated as a rectangular area like the first operation gestalt, but is calculated as a trapezoid area. 
That is, mt- (delta tdwn-f delta tioc) / 2 are used instead of mt-delta tdvm in the formula (6) in the first 
operation gestalt. deltatdwn is the time amoimt dxuing the minimum value time of day tmin when 
taking the maximum time of day tmax and the minimum value Pmmin as mentioned above, in case 
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the pressiire of inhalation of air takes Maximum Pmmax here (deltatdwn=tmin-tmax). deltatioc is the 
time amount between the time of day (valve-opening stage) tio which the cylinder [ i-th ] inlet valve 
6 opens, and the time of day (clausilium stage) tic which an inlet valve 6 closes, i.e., the time amount 
which the cylinder [ i-th ] inlet valve 6 is opening, (deltatioc=tic-tio). Therefore, with the fourth 
operation gestalt, the above-mentioned formula (6) is rewritten like the following formula (12), and 
is used. 
[Equation 10] 



' AD ^ /(Ra Tm\ ^ Atdwn-^Ati 
Mci = tsFmdtwn + mr •■• (12) 

l\Vm) 2 



[0076] 

By the formula (12), the capacity with which the term containing deltaPmdwn is equivalent to the 
area A in drawing 10 (b) is expressed, and since the capacity with which the term containing mt is 
equivalent to the area B in drawing 10 (b) is expressed, the cylinder [ i-th ] restoration air content 
Mci in a cylinder serves as a value adding the area A and area B in drawing 10 (b). The great portion 
of capacity equivalent to the area C shown in drawing 9 can be included in the restoration air content 
in a cylinder by calculating capacity other than the capacity equivalent to area B as a trapezoid so 
that drawing 10 (a) may show. Therefore, according to this operation gestalt, the restoration air 
content Mci in a cylinder serves as a value which expresses more correctly the capacity with which it 
filled up in the cylinder [ i-th ] combustion chamber 5 during the cylinder [ i-th ] valve-opening 
period of an inlet valve 6, and even if it is the case where the valve-opening stage of the inlet valve 6 
of each gas column overlaps, the presumed error of the restoration air content Mci in a cylinder can 
be suppressed small. 

In addition, you may make it the control unit of the fourth operation gestalt calculate the restoration 
air content in an average cylinder combining the control imit of not only the first operation gestalt 
but the second operation gestalt. 
[0077] 

Next, the control unit of the internal combustion engine of the fifth operation gestalt of this invention 
is explained. The control unit of the fifth operation gestalt is the same as the control imit of the first 
operation gestalt fundamentally. By the way, when the operation parameter of internal combustion 
engines, such as an engine rotational frequency, a phase angle of an inlet valve, and a pressure of 
inhalation of air, is the same respectively, the cylinder [ i-th ] restoration air content Mci in a 
cylinder is determined as a meaning with the working angle of the cylinder [ i-th ] inlet valve 6. For 
example, when operation parameters other than a working angle are fixed, the relation between the 
restoration air content Mci in a cylinder and an actual working angle serves as a curve as shown in 
drawing 1 1 . Operation parameters other than a working angle in therefore, the condition of having 
fixed to the specific value or the value of the near, respectively the [ above-mentioned ] — the 
actuation ("air content presumption actuation" is called hereafter) in the third operation gestalt from 
the 1 first operation gestalt — based on a map as shown in drawing 1 1 , the actual working angle of 
the cylinder [ i-th ] inlet valve 6 can be presvimed from the cylinder [ which was therefore presumed / 
i-th ] restoration air content Mci in a cylinder, 
[0078] 

With this operation gestalt, it asks for a curve as shown in drawing 1 1 of the restoration air content in 
a cylinder and a working angle in case operation parameters other than a working angle (for 
example, an engine rotational frequency, the phase angle of an inlet valve 6, the average of the 
pressure of inhalation of air) take a specific value or the value of the near by count experimentally 
beforehand, and, specifically, saves as a map as shown in drawing 1 1 at ROM34 of ECUS 1. And 
when operation parameters other than a working angle have taken the above-mentioned specific 
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value or the value of the near during operation of an internal combustion engine, the restoration air 
content Mci in a cylinder to each gcis column is presumed by the above-mentioned air content 
presumption actuation. The actual working angle of an inlet valve 6 is computed from the restoration 
air content Mci in a cylinder to each presumed gas column, and the map saved at the above ROM 34. 
Thereby, according to this operation gestalt, the actual working angle of an inlet valve 6 can be 
computed comparatively correctly. 
[0079] 

By the way, since an inlet valve 6 is driven, when the good fluctuation valve system (not shown) of 
the electromagnetic inlet valve 6 is prepared, a gap will arise between the target working angle 
ordered by degradation of the spring used into a good fluctuation valve system etc. from ECU31 to a 
good fluctuation valve system, and the actual working angle of an inlet valve 6. Moreover, when the 
inlet valve 6 is driving by the mechanical good fluctuation valve system, a gap will arise between the 
target working angle ordered by wear of the cam used for a good fluctuation valve system etc. from 
ECU31 to a good fluctuation valve system, and the actual working angle of an inlet valve 6. if such a 
gap arises — the working angle of an inlet valve 6 — suitable — being uncontrollable — becoming — 
an engine output, fuel consxmiption, or exhaust air ~ aggravation of description will be caused. 
[0080] 

So, when the actual working angle presumed based on air content presumption actuation differs from 
the target working angle ordered from ECU3 1 to the good fluctuation valve system, it is made to 
always make the actual working angle of an inlet valve 6 in agreement [ by carrying out amendment 
which compensates the difference of the actual working angle and target working angle which were 
presumed ] with a target working angle with this operation gestalt, 
[0081] 

For example, such difference is computed when the actual working angle and target working angle 
which were presumed differ from each other. And a target working angle is ordered from ECU3 1 on 
and after next time in the value to which the difference by which calculation was carried out [ above- 
mentioned ] was added to a good fluctuation valve system. 
[0082] 

therefore, the thing which is controlled to always make the actual working angle of an inlet valve 6 
in agreement with a target working angle according to the fifth operation gestalt — an engine output, 
fuel consumption, or exhaust air ~ aggravation of description can be controlled. 
[0083] 

Next, the control unit of the intemal combustion engine of the sixth operation gestalt of this 
invention is explained. The control xmit of the intemal combustion engine of the sixth operation 
gestalt is the same as that of the first operation gestalt fundamentally. 
[0084] 

By the way, in the third operation gestalt, the restoration air content in a cylinder which can be 
presumed by air content presumption actuation is a thing in front of 1 cycle based on the output from 
the pressure sensor 40 of inhalation of air from the first operation gestalt. That is, with the above- 
mentioned operation gestalt, fuel oil consumption etc. is computed based on the restoration air 
content in a cylinder in front of 1 cycle. Since presumption of the restoration air content in a cylinder 
is performed after filling up with inhalation-of-air gas in a cylinder completely, this is because fuel 
oil consumption etc. cannot be determined based on the restoration air content in a cylinder in the 
same cycle as having presumed the restoration air content in a cylinder. Therefore, from the first 
operation gestalt, in the third operation gestalt, fluctuation of the restoration air content in a cylinder 
is small between cycles, or when there is almost nothing (i.e., only when engine operational status is 
in a steady state), fiiel oil consxunption etc. can be determined based on the computed restoration air 
content in a cylinder. However, when engine operational status is in a transient and it changes the 
restoration air content in a cylinder sharply between cycles, the restoration air content in a cylinder 
presumed by air content presumption actuation in the third operation gestalt from the above- 
mentioned first operation gestalt cannot be used. When engine operational status is in a transient, it is 
necessary to predict the following restoration air content in a cylinder. 
[0085] 

In order to perform such prediction, the restoration air content model Ml 0 in a cylinder mentioned 
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later is used. Although it can predict that the restoration air content in a cylinder of degree cycle 
("the future restoration air content in a cylinder" is called hereafter) mentions later in this restoration 
air content model MIO in a cylinder, the restoration air content in a cylinder computed is not the air 
content for every gas column but the average ("future restoration air content Mcin average cylinder*" 
is called hereafter) of the restoration air content in a cylinder in all gas columns. 
[0086] 

So, with this operation gestalt, fiiture restoration air content Mciin cylinder* for every gas column is 
computed by amending fixture restoration air content Mcin average cylinder' computed with the 
restoration air content model MIO in a cylinder mentioned later. 
[0087] 

The average of all the gas columns of the restoration air content in a cylinder presumed by air 
content presumption actuation in the third operation gestalt from the first operation gestalt is 
specifically computed, and the deflection of the restoration air content in a cylinder of each gas 
column to the average of all these gas columns is computed as correction factor etai. That is, 
correction factor etai to the i-th cylinder is the value which did the division of the restoration air 
content Mci in [ of the i-th cylinder ] a cylinder presumed by the above-mentioned air content 
presumption actuation by the average of all the gas colxmms of the restoration air content in a 
cylinder, as shown in the following formula (13). 
[Equation 11] 



^^-^ ; •••(13) 

^Mci/Ncyl 



In addition, in a formula (13), Ncyl is the number of gas columns. Moreover, sigmaMci is the total 
restoration air content in a cylinder to all the gas columns in 1 cycle, continues for 1 cycle and totals 
the restoration air content Mci in a cylinder presumed by the above-mentioned air content 
presumption actuation. 
[0088] 

And the value which carried out the multiplication of the correction factor etai to the i-th cylinder to 
future restoration air content Mcin average cylinder* computed with the restoration air content model 
MIO in a cylinder mentioned later is made into restoration air content Mciin cylinder [ i-th ] future 
cyUnder* (Mci'=etai-Mc'). Also when future restoration air content Mciin cylinder* can be correctly 
presumed now for every gas column in consideration of dispersion in the restoration air content in a 
cylinder between each gas colvimn and engine operational status is in a transient by this, it can 
perform maintaining the air- fuel ratio of the gaseous mixture in each gas column to a target air-fuel 
ratio etc. In addition, correction factor etai is serially updated, when engine operational status is in a 
steady state, and when engine operational status is in a transient, it is considered as as [ the value 
updated at the end in the last steady state ]. This is because the presumed precision of the restoration 
air content in a cylinder in the transient by the third operation gestalt is low from the first operation 
gestalt. 
[0089] 

In addition, it is good also considering correction factor etai [ as opposed to the i-th cylinder at the 
above-mentioned sixth operation gestalt ] as the average or the weighted average efficiency of a 
correction factor between two or more cycles. For example, weighted-average-efficiency etaiave of a 
correction factor is computed by the following type (14). 
[Equation 12] 
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riiave = s • rji^^^ + (1 - 5) • - (14) 



Here, etai (n) is the correction factor now computed by the formula (13) in the cycle, and etai (n-1) is 
the correction factor computed by the formula (13) in the front cycle. Moreover, s is the weight of a 
weighted average and is a value with which 0<=s<=l is filled and which was defined beforehand. 
Thus, by taking the average or the weighted average efficiency of a correction factor, the error 
produced by the noise of the pressure sensor 40 of inhalation of air etc. can be compensated. 
[0090] 

Drawing 12 shows the operating procediwe which presvunes restoration air content Mciin cylinder 
[ i-th ] fiiture cylinder'. This actuation is performed for every gas colximn. First, in step 161, it is 
judged whether current engine operational status is a steady state. When engine operational status is 
judged to be a steady state, the case where operation parameters, such as an engine rotational 
frequency or an engine load, have fallen within the predetermined range in a fixed period is 
mentioned. When judged with there being no engine operational status in a steady state, steps 162- 
165 are not performed. When judged with engine operational status being in a steady state, it 
progresses to step 162. 
[0091] 

Renewal of correction factor etaiave is performed at steps 162-165 . At step 162, the cylinder [ i-th ] 
restoration air content Mci in a cylinder is presumed by the above-mentioned air content 
presumption actuation. Subsequently, in step 163, the total restoration air content sigmaMci in a 
cylinder to all the gas columns in 1 cycle is computed by adding the cyUnder [ which was computed 
at step 162 / i-th ] restoration air content Mci in a cylinder. Subsequently, at step 164, correction 
factor etai to the i-th cylinder is computed by the formula (13) based on Mci presumed in step 162, 
and sigmaMci computed in step 163. At step 165, weighted-average-efficiency etaiave of a 
correction factor to the i-th cylinder is computed by the formula (14) based on correction factor etai 
computed in step 164 this time and last time. 
[0092] 

Subsequently, at step 166, fiiture restoration air content Mcin average cylinder* computed with the 
restoration air content model Ml 0 in a cylinder is acquired. And in step 167, fiiture restoration air 
content Mcin average cylinder' is made into restoration air content Mciin cylinder [ i-th ] fiiture 
cylinder' by carrying out multiplication at weighted-average-efficiency etaiave of the correction 
factor computed in step 165 (Mci'=etaiave-Mc'). 
[0093] 

Next, the control unit of the intemal combustion engine of the seventh operation gestalt of this 
invention is explained. Although the control unit of the seventh operation gestalt is the same as the 
control xmit of the sixth operation gestalt fimdamentally As opposed to having carried out the 
multiplication of the correction factor etai for every gas column to fiiture restoration air content Mcin 
average cylinder' computed with the restoration air content model Ml 0 in a cylinder later mentioned 
with the sixth operation gestalt With this operation gestalt, fiiture restoration air content Mciin 
cylinder' for every gas column is computed by adding amendment capacity deltaMi for every gas 
column to fiiture above-mentioned restoration air content Mcin average 
cylinder' (Mci'=Mc'+deltaMi). 
[0094] 

Here, the calculation approach of amendment capacity deltaMi is explained. The deflection of the 
restoration air content in a cylinder of each gas column to the restoration air content in an average 
cylinder changes according to the value of an intemal combustion engine's operation parameter (for 
example, a working angle, an engine rotational frequency, a phase angle). For example, when a 
working angle is large about the same gas column, deflection is small, and when a working angle is 



http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_web_cgi__ejje 9/2 1/2006 



IP,2005-090437,A [DETAILED DESCRIPTION] 



Page 19 of 27 



small, deflection is large, when a working angle is taken for an example. Since amendment capacity 
deltaMi is for compensating this deflection, it needs to determine that it becomes this deflection and 
the same value. Therefore, when the working angle VL is large, it is necessary to set up amendment 
capacity, and as shown in drawing 13 , when the working angle VL is small, it is necessary, as for 
the relation between the working angle VL and amendment capacity deltaMi, to set up amendment 
capacity greatly small. 
[0095] 

Moreover, the relation between the working angle VL and the above-mentioned deflection changes 
with extent, such as between gas colxmms and aging. Therefore, the relation between the working 
angle VL and amendment capacity deltaMi turns into various relation between gas columns and by 
aging etc. similarly like o, p, and q which were shown in drawing 13 . 
[0096] 

So, with this operation gestalt, first, it asks for the relation between the working angle VL and 
amendment capacity deltaMi experimentally beforehand, and saves as a map at ROM34 of ECU31. 
And the working angle VL and the cylinder [ i-th ] restoration air content Mci in a cylinder at the 
time at the time of some detection conditions are presumed by air content presumption actuation of 
the first operation gestalt to the third operation gestalt. and future restoration air content Mcin 
average cylinder' computed with the restoration air content model Ml 0 in a cylinder in this cycle 
from the cylinder [ which was presimied / i-th ] restoration air content Mci in a cylinder — 
subtracting — the account of a top — amendment capacity deltaMi in some detection conditions is 
computed. For example, when the amendment capacity computed when a working angle was VLl is 
deltaMi 1, as shown in drawing 13 , this point is on Curve o, therefore Curve o is adopted as a ciuve 
of the amendment capacity about the i-th cylinder. 
[0097] 

And by the i-th cylinder, amendment capacity deltaMi is computed from the map shown at drawing 
13 based on the working angle VL for every cycle after degree cycle, and let fixture restoration air 
content Mciin cylinder' the i-th cylinder be the value with which the above-mentioned amendment 
capacity deltaMi was added to restoration air content Mcin average cylinder' (Mci'=Mc'+deltaMi). 
[0098] 

Such actuation is performed for every gas column, thereby, dispersion in the restoration air content 
in a cylinder between each gas column can be compensated, and future restoration air content Mciin 
cylinder' can be correctly computed for every gas column. 
[0099] 

Next, the restoration air content model Ml 0 in a cylinder is explained. In addition, below, Mc' and 
the intake air flow in an average cylinder are made into mc' for the restoration air content in an 
average cylinder computed with the restoration air content model Ml 0 in a cylinder. 
[0100] 

The restoration air content model Ml 0 in a cylinder is equipped with the electronics control throttle 
model Ml 1, the throttle model M12, the inlet-pipe model 1VI13, and the inlet-valve model M14 as 
shown in drawing 14 . The accelerator pedal control input Accp detected by the load sensor 46 is 
inputted into the electronics control throttle model Ml 1, and throttle opening ("read-ahead throttle 
opening" is called hereafter) thetat to which the actual throttle valve 18 reaches after predetermined 
time deltaT is outputted to it. Read-ahead throttle opening thetat outputted to the throttle model Ml 2 
fi-om the electronics control throttle model Mil, Atmospheric pressure of the perimeter of an 
internal combustion engine detected by the atmospheric pressure sensor 44 0 Or the pressure Pa of 
the air inhaled by the inlet pipe 15 and atmospheric temperature of the perimeter of an internal 
combustion engine detected by the atmospheric temperature sensor 43 () Or the temperature Ta of 
the air inhaled by the inlet pipe 15 and the pressure Pm in the inhalation-of-air branch pipe 13 
computed in the inlet-pipe model Ml 3 mentioned later (pressure of inhalation of air) are inputted. By 
substituting for the model type of the throttle model M12 which mentions the value of each these- 
inputted parameter later, the flow rate ("the throttle- valve passage air flow rate mt" is called 
hereafter) of the air which passes a throttle valve 18 to per unit time amount is computed. The 
throttle-valve passage air flow rate mt computed in the throttle model Ml 2 is inputted into the inlet- 
pipe model Ml 3. 
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[0101] 

The throttle- valve passage air flow rate mt computed in the throttle model Ml 2 by the inlet-pipe 
model Ml 3, and the flow rate of the inhalation-of-air gas which flows into per [ which is explained 
in full detail below ] unit time amount in a combustion chamber 5 ("intake-air-flow mcin average 
cylinder'" is called hereafter.) in addition, about the definition of intake-air-flow mcin average 
cylinder' the inlet- valve model Ml 4 — setting — explaining in full detail — by substituting for the 
model type of the inlet-pipe model Ml 3 which is inputted and mentions the value of each these- 
inputted parameter later The temperature ("the temperature Tm of inhalation of air" is called 
hereafter) of the inhalation-of-air gas which exists in the pressure ("the pressure Pm of inhalation of 
air" is called hereafter) of the inhalation-of-air gas which exists in the inhalation-of-air branch pipe 
13 and a surge tank 14, the inhalation-of-air branch pipe 13, and a surge tank 14 is computed. Both 
the pressures Pm of inhalation of air and temperature Tm of inhalation of air that were computed in 
the inlet-pipe model Ml 3 are inputted into the inlet- valve model Ml 4, and the pressxire Pm of 
inhalation of air is further inputted also into the throttle model Ml 2. 
[0102] 

Intake-air-flow mcin average cylinder* is computed by the inlet-valve model M14 by substituting for 
the model type of the inlet- valve model Ml 4 with which the atmospheric temperature Ta other than 
the pressure Pm of inhalation of air computed in the inlet-pipe model Ml 3 and the temperature Tm 
of inhalation of air is inputted, and mentions the value of each these-inputted parameter later. 
Computed intake-air- flow mcin average cylinder' is changed into restoration air content Mcin 
average cylinder', and the fuel oil consumption firom a fuel injection valve is determined based on 
this restoration air content Mcin average cylinder'. Moreover, intake-air-flow mcin average cylinder' 
computed in the inlet-pipe model Ml 3 is inputted into the inlet-pipe model Ml 3, 
[0103] 

Since the value of the parameter computed in a certain model in the restoration air content model 
MIO in a cylinder is used as an input value to another model so that dravsdng 14 may show, in the 
restoration air content model MIO in a cylinder whole, the values actually inputted are only three 
parameters, throttle opening thetat, atmospheric pressure Pa, and atmospheric temperature Ta, and 
restoration air content Mcin average cylinder' is computed from these three parameters. 
[0104] 

Next, each models Ml 1 -Ml 4 of the restoration air content model MIO in a cylinder are explained. 

[0105] 

The electronics control throttle model Mil is a model which presumes throttle opening ("read-ahead 
throttle opening" is called hereafter) thetat to which the actual throttle valve 18 reaches after 
predetermined time deltaT based on the accelerator pedal control input Accp detected by the load 
sensor 46. In this operation gestalt, throttle opening thetat is called for based on the map which 
specifies the relation between the accelerator pedal control input Accp shown in the accelerator pedal 
control input Accp detected by the load sensor 46 by throttle- valve electronics control logic, and 
drawing 15 , and target throttle opening thetat. Thus, called- for throttle opening thetat is sent to the 
throttle model Ml 2. On the other hand, it asks as target throttle opening thetar with the final value 
fi-om which only predetermined time deltaT (for example, 64msec) delayed this throttle opening 
thetat, and a driving signal is sent out to a step motor 17 so that the actual throttle opening TA may 
be set to target throttle opening tiietar. 
[0106] 

Thus, target throttle opening thetar is equal to throttle opening thetat as which only predetermined 
time deltaT was determined according to the accelerator pedal control input Accp at the last time 
firom this time. Since a throttle valve 1 8 drives based on target throttle opening thetar, as for throttle 
opening thetat, only deltaT is previous throttle opening from the throttle opening of the actual 
throttle valve 18. Conversely, if it says, throttle opening thetat is the throttle opening to which the 
actual throttle valve 18 reaches after predetermined time deltaT. 
[0107] 

In the throttle model Ml 2, the throttle- valve passage air flow rate mt is computed based on the 
following formula (15) from read-ahead throttle opening thetat outputted from an atmospheric 
pressure Pa, atmospheric temperature Ta, the pressure Pm of inhalation of air, and the electronics 
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control throttle model Mil. Here, mu in a formula (15) is a flow coefiBcient in a throttle valve, is the 
function of throttle opening thetat and becomes settled from a map as therefore shown in drawing 
16 . Moreover, At shows the opening cross section of a throttle valve, is the function of throttle 
opening thetat and becomes settled from a map as shown in drawing 17 . In addition, mu-At which 
summarized these flow coefficients mu and the opening cross section At may be made to be found 
on one map from throttle opening thetat. Moreover, Ra is a constant about a gas constant and is the 
value which did the division of the gas constant with the mass Mlmol of the gas per mol (air) in fact. 
[Equation 13] 



mt^fi^At' ^ — _ >ol-^j •••(15) 



Pa 
^RaTa 



[0108] 

Moreover, phi (Pm/Pa) is the function shown in the following type (16), and kappa in this formula 
(16) is the ratio of specific beat (it considers as constant value). Since this function phi (Pm/Pa) can 
be expressed with a graph as shown in drawing 18 , it saves by using such a graph as a map at 
ROM34 of ECU31, and does not calculate using a formula (16) in fact, but you may make it 
calculate the value of phi (Pm/Pa) from a map. 
[Equation 14] 



2(it4.l) ^'"^ 1 




Pa ) Pa\\ Pa) ^^ + 1 



Pm\^Pm 



Pa ^ —(16) 



[0109] 

The formula (15) and formula (16) of these throttle model M12 Atmospheric temperature Ta and the 
pressure of the gas of the lower stream of a river of a throttle valve 1 8 are made [ the pressure of the 
gas of the throttle- valve 1 8 upstream ] into the pressure Pm of inhalation of air for the temperature of 
the gas of atmospheric pressure Pa and the throttle- valve 18 upstream. To the model of a throttle 
valve 18 as shown in drawing 19 , the conservation of mass, the law of conservation of energy, and a 
law of conservation of momentum are applied, and it is obtained by using a gaseous equation of 
state, the definition equation of the ratio of specific beat, and a Mayer's relation equation further. 
[0110] 

In the inlet-pipe model M13, the pressure Pm of inhalation of air and the temperature Tm of 
inhalation of air are computed based on the following formula (17) and a formula (18) from the 
throttle-valve passage air flow rate mt, intake-air- flow mcin average cylinder*, and atmospheric 
temperature Ta. In addition, Vm in a formula (17) and a formula (18) is a constant equal to the 
volume of the part ("a part for an inhalation-of-air tube part" is called hereafter) of the inlet-pipe 13 
grade from a throttle valve 18 to an inlet valve 6. 
[Equation 15] 
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d(Pm\ R , 



dPm Ra , ^ ^ ^ 

—— = K (mtTa-mcTm) -(18) 

dt Vm 



[0111] 

Here, the inlet-pipe model Ml 3 is explained with reference to drawing 20 . When the total quantity 
of gas for an inhalation-of-air tube part (the total inhalation-of-air capacity) is set to M, the temporal 
response of the total quantity of gas M Since it is equal to the flow rate mt of the gas which flows 
into a part for an inhalation-of-air tube part, i.e., a throttle- valve passage air flow rate, and the flow 
rate of the gas which flows out of a part for an inhalation-of-air tube part, i.e., a difference with 
intake-air-flow mcin average cylinder'. The following equation (19) is obtained by the conservation 
of mass, and an equation (17) is obtained from this equation (19) and a gaseous equation of state 
(Pm-Vm=M-R-Tm). 
[Equation 16] 



^mt-mc -"(19) 

dt 



[0112] 

Moreover, the amount of temporal responses of energy M-Cv-Tm of the gas for an inhalation-of-air 
tube part is equal to the difference of the energy of the gas which flows into a part for an inhalation- 
of-air tube part, and the energy of the gas which flows out of a part for an inhalation-of-air tube part. 
For this reason, if temperature of the gas which flows out of a part for atmospheric temperature Ta 
and an inhalation-of-air tube part the temperature of the gas which flows into a part for an inhalation- 
of-air tube part is made into the temperature Tm of inhalation of air, the following equation (20) will 
be obtained by the law of conservation of energy, and an equation (18) will be obtained from this 
equation (20) and the equation of state of the above-mentioned gas. 
[Equation 17] 



^^^^^""'^^^ =Cp^mt'Ta-CpmC'Tm -(20) 
dt 



[0113] 

In the inlet- valve model Ml 4, intake-air-flow mcin average cylinder' is computed based on the 
following type (21) from the pressure Pm of inhalation of air, the temperature Tm of inhalation of 
air, and atmospheric temperature Ta. In addition, a in a formula (21) and b are values which become 
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settled from the phase angle of an inlet valve 6, and a working angle in the case of the internal 
combustion engine which had the good fluctuation valve system which can change the phase angle 
(valve timing) and working angle of an inlet valve 6 further from the engine engine speed Ne. 
[Equation 18] 



mc = • (a ' Pm - b) 



-(21) 



[0114] 

The inlet-valve model Ml 4 mentioned above is explained with reference to drav^ng 21 . It decides, 
when an inlet valve 6 closes the valve (at the time of an inhalation-of-air valve-closing valve), and 
restoration air content Mcin average cylinder* which is the amoimt of the inhalation-of-air gas 
currently generally inhaled in the combustion chamber 5 when an inlet valve 6 closes is proportional 
to the pressure in the combustion chamber 5 at the time of an inhalation-of-air valve-closing valve. 
Moreover, it can be considered that the pressure in the combustion chamber 5 at the time of an 
inhalation-of-air valve-closing valve is equal to the pressure Pm of the gas of the inlet- valve 
upstream, i.e., the pressure of inhalation of air. Therefore, restoration air content Mcin average 
cylinder* can be approximated if it is proportional to the pressure Pm of inhalation of air. 
[0115] 

What equalized the amount of all the inhalation-of-air gas that flows into per unit time amount out of 
a part for an inhalation-of-air tube part here, Or if what continued like the inhalation-of-air line of 
one g£is column, and equalized the amount of the inhalation-of-air gas inhaled in all the combustion 
chambers 5 from a part for an inhalation-of-air tube part (a part for 180 degrees of crank angles [ It 
mentions later, like this operation gestalt ]) is made into intake-air-flow mcin average cylinder* (it 
explains in full detail below) per unit time amoxmt Since restoration air content Mcin average 
cylinder' is proportional to the pressure Pm of inhalation of air, it is thought that intake-air-flow mcin 
average cylinder' is also proportional to the pressure Pm of inhalation of air. Based on the theory and 
a mle of thumb, the above-mentioned formula (21) is obtained from this. In addition, the value a in a 
formula (21) is a proportionality coefficient, and can be foimd from the three-dimensions map which 
made the parameter the engine rotational frequency Ne, the amount indicated value VL of lifts of an 
inlet valve 6, and phase angle indicated value VT of an inlet valve 6. In addition, this three- 
dimensions map is beforehand called for by coimt experimentally, and is saved at ROM34 of 
ECUS 1 . A value b is a value (it is thought that the division was done at deltaTl 80 degrees of time 
amount which mentions later the amount of burnt gases which remains in a combustion chamber 5 at 
the time of exhaust valve 8 clausilium) showing the burnt gas which remained in the combustion 
chamber 5. Moreover, in actual operation, since the temperature Tm of inhalation of air may change 
to a transient a lot, the multiplication of Ta/Tm drawn based on the theory and a rule of thumb as 
amendment to this is carried out. 
[0116] 

Here, with reference to [* / intake-air- flow mcin average cylinder ] drawing 22 , the case where an 
internal combustion engine is a 4-cylinder is explained. In addition, drawing 22 is the flow rate of 
the inhalation-of-air gas by which an axis of abscissa flows into angle of rotation of a crankshaft per 
unit time amount, and an axis of ordinate actually flows into a combustion chamber 5 from a part for 
an inhalation-of-air tube part. As shown in drawing 22 , in the internal combustion engine of a 4- 
cylinder, it opens in order of a No. 1 gas column, a No. 3 gas column, a No. 4 gas colimin, and a No. 
2 gas column, and inhalation-of-air gas flows [ an inlet valve 6 ] into the combustion chamber 5 of 
each gas column from a part for an inhalation-of-air tube part according to the amoimt of valve 
opening of the inlet valve 6 corresponding to each gas column. For example, the variation rate of the 
flow rate of the inhalation-of-air gas which flows in the combustion chamber 5 of each gas column 



http://www4.ipdl.ncipi.go.jp/cgi-bin/tran__web_cgi_ejje 



9/21/2006 



IP,2005-090437,A [DETAILED DESCRIPTION] 



Page 24 of 27 



from a part for an inhalation-of-air tube part is as the broken line having shown to drawing 22 , and 
the flow rate of the inhalation-of-air gas which synthesizes this and flows into the combustion 
chamber of all gas columns from the inhalation-of-air branch pipe 13 is as the continuous line having 
shown to drawing 22 . Moreover, restoration air content Mcin average cylinder* to a No. 1 gas 
column is as the slash having shown to drawing 22 , for example. 
[0117] 

On the other hand, it is intake-air-flow mcin average cylinder* which equalized the amoimt of the 
inhalation-of-air gas which flows into the combiistion chemiber of all gas columns from the inlet pipe 
shown as the continuous line, and an altemate long and short dash line shows all over drawing. And 
in the case of a 4-cylinder, what carried out the multiplication of the deltaTlSO degrees of the time 
amount concerning a crankshaft carrying out 180-degree (namely, include angle which divided 
include angle of 720 degrees which crankshaft rotates in 1 cycle in 4 stroke type internal combustion 
engine by nxmiber of gas columns) rotation becomes restoration air content Mcin average cylinder' at 
intake-air-flow mcin average cylinder* shown with this altemate long and short dash line. TTierefore, 
restoration air content Mcin average cylinder' is computed by carrying out the multiplication of the 
deltaTlSO degrees to intake-air-flow mcin average cylinder' computed with the inlet-valve model 
Ml 4 (Mc'=mc', deltaTlSO degrees). In consideration of restoration air content Mcin average 
cylinder* being proportional to the pressure at the time of an inhalation-of-air valve-closing valve, 
what carried out the multiplication of the deltaTlSO degrees to intake-air-flow mcin average 
cylinder* at the time of an inhalation-of-air valve-closing valve is made more into restoration air 
content Mcin average cylinder* at a detail. 
[0118] 

Next, the above-mentioned restoration air content model Ml 0 in a cylinder is mounted in an internal 
combustion engine's control imit, and the case where restoration air content Mcin average cylinder* is 
actually computed is explained. Restoration air content Mcin average cylinder' is expressed using the 
restoration air content model 10 in a cylinder by solving the above-mentioned formula (15), a 
formula (17), a formula (18), and a formula (21). In this case, these formulas need to be discretized 
in order to process by ECU31. If a formula (15), a formula (17), a formula (IS), and a formula (21) 
are discretized using time of day t and count spacing deltat, the following type (22), a formula (23), a 
formula (24), and a formula (25) will be obtained, respectively. In addition, the temperature Tm 
(t+deltat) of inhalation of air is computed by the formula (26) from Pm/Tm (t+deltat) computed by 
the formula (23) and the formula (24), respectively and Pm (t+deltat). 
[Equation 19] . 
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^R Ta \ Pa ) 



(22) 



— (/ + AO = —(0 + A/ -(m/CO - mc{t)) - (23) 

Tm Tm Vm 

n 

Pm{t + A/) = Pm{t) + A/ • ^ imt{t) - Ta - mc(0 ■ Tm{t)) — (24) 

mc(0 = • (a • />m(0 --b) - (25) 

rm(f) 

7)«(> + A/)= /"C^^) ...(26) 



[0119] 

thus, in the mounted restoration air content model Ml 0 in a cylinder The throttle- valve passage air 
flow rate mt in the time of day t computed by the formula (22) of the throttle model Ml 2 (t), Litake- 
air-flow mcin average cylinder* (t) in the time of day t computed by the formula (25) of the inlet- 
valve model Ml 4 It is substituted for the formula (23) and formula (24) of the inlet-pipe model Ml 3, 
and, thereby, the pressure Pm of inhalation of air in time-of-day t+delta t (t+deltat) and the 
temperature Tm (t+deltat) of inhalation of air are computed. Subsequently, Pm (t+deltat) and Tm 
(t+deltat) which were computed are substituted for the formula (22) and formula (25) of the throttle 
model Ml 2 and the inlet- valve model Ml 4, and, thereby, the throttle- valve passage air flow rate mt 
(t+deltat) in time-of-day t+delta t arid iiitalce-air-flow mcin average cylinder' (t+deltat) are computed, 
and ~ being such — coimt ~ repeating — things ~ a read ahead — a throttle — opening — theta ~ t ~ 
atmospheric pressure — Pa — and — atmospheric temperature ~ Ta — from ~ arbitration ~ time of 
day — t — it can set — an average — a cylinder — inside — an intake air flow — mc — ' — computing ~ 
having — computing — having had — an average — a cylinder — inside — an intake air flow — mc — ' - 
- the above — time amoimt ~ delta — T ~ 180 ~ degree — multiplication ~ carrying out ~ things — 
arbitration — time of day — t — it can set — an average — a cylinder — inside -- restoration — an air 
content — Mc — * — computing — having . Since read-ahead throttle opening thetat is throttle opening 
only with early deltaT from the throttle opening of the actual throttle valve 18 especially, restoration 
air content Mcin average cylinder' computed is also a ftiture value. 
[0120] 

In addition, in the time t= 0 of an internal combustion engine's starting, i.e., time of day, it is 
supposed that the pressure Pm of inhalation of air is equal to atmospheric pressure (Pm(0) =Pa), it is 
supposed that the temperature Tm of inhalation of air is equal to atmospheric temperature (Tm(0) 
=Ta), and the count in each models Ml 1-M13 is started. 
[0121] 

In addition, although [ the above-mentioned restoration air content model Ml 0 in a cylinder ] 
atmospheric temperature Ta and atmospheric pressure Pa are fixed The value detected in time of day 
t by the atmospheric temperature sensor for it being good also as a value which changes with time of 
day, for example, detecting atmospheric temperature Atmospheric temperature Ta (t) You may make 
it substitute for the above-mentioned formula (22), a formula (24), and a formula (25) by making 

http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_web_cgi_ejje 9/2 1/2006 



JP,2005.090437,A pETAILED DESCRIPTION] 



Page 26 of 27 



into atmospheric pressure Pa (t) the value detected in time of day t by the atmospheric pressure 

sensor for detecting atmospheric pressure. 

[0122] 

In addition, in this specification, the case where engine operational status is in a steady state means 
operational status which is maintained almost uniformly, without hardly changing an intemal 
combustion engine's operation parameters (for example, an engine rotational frequency, an engine 
load, the restoration air content in a cyUnder, etc.), and, on the other hand, the case where engine 
operational status is in a transient means operational status to which an intemal combustion engine's 
operation parameter is changed sharply 
[Brief Description of the Drawings] 
[0123] 

[Drawing 1] It is drawing showing the whole intemal combustion engine by which the control unit of 
the intemal combustion engine of this invention is used. 

[Dravying 2] It is drawing showing the fundamental concept of the inlet-pipe model of this invention. 

[Drawing 3] It is drawing showing change of the flow rate to a crank angle, and change of the 
pressure of inhalation of air. 

[Drawing 4] It is the flow chart which shows the operating procedure of presumption of the 
restoration air content in a cylinder to each gas column. 

[Drawing 5] The relation of the ignition timing and combustion energy in each gas column is shown. 

[Drawdng 6] It is the flow chart which shows the operating procedure which sets the fiiel oil 
consumption and ignition timing in each gas column. 

[Drav^ing 7] It is the flow chart which shows the operating procedure which equalizes and presumes 
the restoration air content in a cylinder to each gas column between cycles. 
[Drawing 8] It is the same drawing as drawing 3 for explaining the presumed approach of the 
restoration air content in a cylinder in the third operation gestalt. 

[Drawing 9] It is drawing showing change of the flow rate to the crank angle in the case of 
overlapping the valve-opening stage of the inlet valve 6 between gas columns. 
[Drawing 10] It is drawing showing change of the flow rate to a crank angle. 

[Dravying 11] It is drawing showing the relation between the amount of filler gas in a cylinder, and a 
working angle. 

[Drawing 12] It is the flow chart which shows the operating procedure which presvmies restoration 
air content Mciin cylinder [ i-th ] future cylinder*. 

[Drawing 13] It is drawing showing the relation between a working angle and amendment capacity 
deltaMi. 

[Drawing 14] It is drawing showing the amount model of inspired gas used by this invention. 
[Drawing 15] It is drawing showing the relation between an accelerator pedal control input and 
target throttle opening. 

[Drawing 16] It is drawing showing the relation between throttle opening and a flow coefficient. 
[Drawing 17] It is drawing showing the relation between throttle opening and the opening cross 
section. 

[Drawing 18] It is drawing showing Function phi (Pm/Pa). 

[Drawing 19] It is drawing showing the fundamental concept of a throttle model. 

[Drawing 20] It is drawing showing the fimdamental concept of an inlet-pipe model. 

[DraMang 21] It is drawing showing the fundamental concept of an inlet- valve model. 

[Drawing 22] It is drawing about the definition of the restoration air content in a cylinder, and the 

inhalation-of-air capacity in a cylinder. 

[Description of Notations] 

[0124] 

1 — Engine body 

5 — Combustion chamber 

6 — Inlet valve 

7 ~ Suction port 
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8 — Exhaust valve 

1 1 — Fuel injection valve 

13 ~ Downstream inlet pipe 

14 -- Surge tank 

1 5 — Upstream inlet pipe 
18 -Throttle valve 

19 — Air flow meter 

31 — ECU (electronic control unit) 

40 — Pressure sensor of inhalation of air 

41 — Temperature sensor of inhalation of air 



[Translation done.] 
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* NOTICES * 
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DESCRIPTION OF DRAWINGS 



[Brief Description of the Drawings] 

[0123] 

[Drawing 1] It is drawing showing the whole intemal combustion engine by which the control unit of 
the intemal combustion engine of this invention is used. 

[Drawing 2] It is drawing showing the fundamental concept of the inlet-pipe model of this invention. 

[Drawing 3] It is drawing showing change of the flow rate to a crank angle, and change of the 
pressure of inhalation of air. 

[Drawing 4] It is the flow chart which shows the operating procedure of presumption of the 
restoration air content in a cylinder to each gas column. 

[Drawing 5] The relation of the ignition timing and combustion energy in each gas column is shown. 

[Drawing 6] It is the flow chart which shows the operating procedure which sets the fuel oil 
consumption and ignition timing in each gas coliram. 

[Drawing 7] It is the flow chart which shows the operating procedure which equalizes and presumes 
the restoration air content in a cylinder to each gas column between cycles. 
[Drawing 8] It is the same drawing as drawing 3 for explaining the presumed approach of the 
restoration air content in a cylinder in the third operation gestalt. 

[Drawing 9] It is drawing showing change of the flow rate to the crank angle in the case of 
overlapping the valve-opening stage of the inlet valve 6 between gas columns. 
[Drawing 10] It is drawing showing change of the flow rate to a crank angle. 

[Drawing 11] It is drawing showing the relation between the amoimt of filler gas in a cyUnder, and a 
working angle. 

[Drawing 12] It is the flow chart which shows the operating procedure which presumes restoration 
air content Mciin cylinder [ i-th ] future cylinder'. 

[Drawing 13] It is drawing showing the relation between a working angle and amendment capacity 
deltaMi. 

[Drawing 14] It is drawing showing the amount model of inspired gas used by this invention. 
[Drawing 15] It is drawing showing the relation between an accelerator pedal control input and 
target throttle opening. 

[Drawing 16] It is drawing showing the relation between throttle opening and a flow coefficient. 
[Drawing 17] It is drawing showing the relation between throttle opening and the opening cross 
section. 

[Drawing 18] It is drawing showing Function phi (Pm/Pa). 

[Drawing 19] It is drawing showing the fundamental concept of a throttle model. 

[Drawing 20] It is drawing showing the fundamental concept of an inlet-pipe model. 

[Drawing 21] It is drawing showing the fundamental concept of an inlet-valve model. 

[Drawing 22] It is drawing about the definition of the restoration air content in a cylinder, and the 

inhalation-of-air capacity in a cylinder. 
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[Drawing 3] 
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[Procedure revision] 





[Filing Date] July 1 1 , Heisei 1 8 (2006. 7.11) 
[Procedure amendment 1] 
[Document to be Amended] Claim 
[Item(s) to be Amended] Claim 1 
[Method of Amendment] Modification 
[The contents of amendment] 
[Claim 1] 

A throttle passage air content calculation means to compute the throttle passage air content which 
passes a throttle valve, 

An excessive air content calculation means to compute the excessive air content to the above- 
mentioned gas column equivalent to the amount of pressure drawdowns accompanjdng pulsation of 
the pressure of inhalation of air by the inlet valve corresponding to each gas column opening, 
A restoration air content presumption means in a cylinder to presume the restoration air content in a 
cylinder for every gas column based on the throttle passage air content detected by the above- 
mentioned throttle passage air content detection means, and the excessive air content computed by 
the above-mentioned excessive air content calculation means. 

The control unit of the internal combustion engine possessing flie engine control means which 
controls an internal combustion engine based on the restoration air content in a cylinder for every gas 
column presxuned by this restoration air content presimiption means in a cylinder. 
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[Procedure amendment 2] 

[Document to be Amended] Specification 

[Item(s) to be Amended] 0008 

[Method of Amendment] Modification 

[The contents of amendment] 

[0008] 

A throttle passage air content calculation means to compute the throttle passage air content which 
passes a throttle valve in the 1st invention in order to solve the above-mentioned technical problem, 
An excessive air content calculation means to compute the excessive air content to the above- 
mentioned gas column equivalent to the amount of pressure drawdowns accompanying pulsation of 
the pressure of inhalation of air by the inlet valve corresponding to each gas column opening, A 
restoration air content presumption means in a cylinder to presume the restoration air content in a 
cylinder for every gas column based on the throttle passage air content detected by the above- 
mentioned throttle passage air content detection means, and the excessive air content computed by 
the above-mentioned excessive air content calculation means, The control unit of the internal 
combustion engine possessing the engine control means which controls an internal combustion 
engine based on the restoration air content in a cylinder for every gas colxmm presumed by this 
restoration air content presumption means in a cylinder is offered. 
[Procedvu'e amendment 3] 
[Document to be Amended] Specification 
[Item(s) to be Amended] 0015 
[Method of Amendment] Modification 
[The contents of amendment] 
[0015] . 

The flow rate sensor which detects the throttle-valve passage air flow rate which passes a throttle 

valve in any 6th one invention is provided, the 8th invention — the 1- The period during the 

minimxmi value stage when the maximvun stage when the pressure of inhalation of air serves as max 

in the period which the inlet valve corresponding to each gas column is opening, and the period of 

the near, and the pressure of inhalation of air serve as min is set to deltatdwn. If the period between 

the valve-opening stage of the above-mentioned inlet valve and a clausilium stage is set to deltatioc 

and the throttle- valve passage air flow rate detected by the above-mentioned flow rate sensor during 

these periods is set to mt The above-mentioned throttle passage air content calculation means 

computes the throttle passage air content Mt based on the following type (1). 

Mt=mt- (delta tdwn+delta tioc)/2 - (1) 

[Procedure amendment 4] 

[Document to be Amended] Specification 

[Item(s) to be Amended] 0033 

[Method of Amendment] Modification 

[The contents of amendment] 

[0033] 

If the inflow air flow rate (mt) for an inhalation-of-air tube part is larger than an outflow air flow rate 
(mci), the pressure of inhalation of air will rise, if small, the pressure of inhalation of air will 
descend, if equal, the pressure of inhalation of air is fixed, and a formula (3) shows that variation 
deltaPm of the pressure of inhalation of air for [ delta ] t seconds is equivalent to a changed part of 
the air content in an inhalation-of-air tube part part. In addition, when it is in a steady state which 
engine operational status mentions later, the amount of inhalation-of-air tube part is buffered to 
being intermittent depending on closing motion of an inlet valve 6, and the change is [ the inflow air 
flow rate (mt) fi-om a part for an inhalation-of-air tube part ] quiet [ the outflow air flow rate (mci) 
fi-om a part for an inhalation-of-air tube part ]. For this reason, the size relation between an outflow 
air flow rate (mci) and an inflow air flow rate (mt) repeats reversal (refer to drawing 3 (a)). This 
means that the value in the parenthesis of the right-hand side of the above-mentioned formula (2) 
repeats positive/negative reversal a fixed period, i.e., the pressure of inhalation of air repeats a rise 
and descent a fixed period, and expresses pulsation of the pressxire of inhalation of air. 
[Procedure amendment 5] 
[Document to be Amended] Specification 



http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_web_cgi_ejje?u=http%3A%2F%2Fwww4,ip... 9/21/2006 



Page 3 of 4 



[Item(s) to be Amended] 0042 
[Method of Amendment] Modification 
[The contents of amendment] 
[0042] 

Subsequently, at step 104, it is judged whether it is in the condition that the value of the closing 

motion valve flag Vlv (n-1) in time of day n-1 is set to 1, and the value of the closing motion valve 

flag Vlv in time of day n (n) is set to 0. That is, at step 104, it is judged whether whether the inlet 

valve's 6 opening at the time of the last actuation, and the inlet valve's 6 closing the valve at the time 

of this actuation and the time of this actuation that is, are [ inlet valves 6 ] the times of clausilium 

being carried out. In step 104, the time of this actuation progresses to step 105, when it judges that an 

inlet valve 6 is not a time of clausiliimi being carried out. 

[Procedure amendment 6] 

[Document to be Amended] Specification 

[Item(s) to be Amended] 0064 

[Method of Amendment] Modification 

[The contents of amendment] 

[0064] 

At step 156, 1 is added to the value eye of a cycle counter. A cycle counter is a counter showing the 
niunber of cycles firom an engine start up. Subsequently, at step 157, the restoration air content Mci 
(eye) in a cylinder in this cycle eye is computed by tiie above-mentioned formula (6) like step 1 16 of 
drawing 4. 

[Procedure amendment 7] 

[Document to be Amended] Specification 

[Item(s) to be Amended] 0075 

[Method of Amendment] Modification 

[The contents of amendment] 

[0075] 

So, with the fourth operation gestalt, it is supposed that any capacity other than the capacity which is 
equivalent to area A among the restoration air contents Mci in a cylinder to each gas column is not 
calculated as a rectangular area like the first operation gestalt, but is calculated as a trapezoid area. 
That is, mt- (delta tdwn+delta tioc) / 2 are used instead of mt-delta tdwn in the formula (6) in the fu-st 
operation gestalt. Here, deltatdwn is the time amount during the minimum value time of day tmin 
when taking the maximum time of day tmax and the minimum value Pmmin as mentioned above, in 
case the pressure of inhalation of air takes Maximum Pmmax (deltatdwn=tmin-tmax), deltatioc is the 
time amount between the time of day (valve-opening stage) tio which the cylinder [ i-th ] inlet valve 
6 opens, and the time of day (clausilium stage) tic which an inlet valve 6 closes, i.e., the time amount 
which the cylinder [ i-th ] inlet valve 6 is opening, (deltatioc=tic-tio). Therefore, with the fourth 
operation gestalt, the above-mentioned formula (6) is rewritten like the following formula (12), and 
is used. 
[Equation 10] 



Mci = IsPrndwrii 



Atdwn + Atioc 
+ mt 



[Procedure amendment 8] 
[Docxmient to be Amended] Specification 
[Item(s) to be Amended] 0076 
[Method of Amendment] Modification 
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[The contents of amendment] 
[0076] 

By the formula (12), the capacity with which the term containing deltaPmdwn is equivalent to the 
area A in drawing 10 (b) is expressed, and since the capacity with which the term containing mt is 
equivalent to the area B in drawing 10 (b) is expressed, the cylinder [ i-th ] restoration air content 
Mci in a cylinder serves as a value adding the area A and area B in drawing 10 (b). The great portion 
of capacity equivalent to the area C shown in drawing 9 can be included in the restoration ak content 
in a cylinder by calculating capacity other than the capacity equivalent to area A as a trapezoid so 
that drawing 10 (a) may show. Therefore, according to this operation gestalt, the restoration air 
content Mci in a cylinder serves as a value which expresses more correctly the capacity with which it 
filled up in the cylinder [ i-th ] combustion chamber 5 during the cyUnder [ i-th ] valve-opening 
period of an inlet valve 6, and even if it is the case where the valve-opening stage of the inlet valve 6 
of each gas coliunn overlaps, the presumed error of the restoration air content Mci in a cylinder can 
be suppressed small. 

In addition, you may make it the control imit of the fourth operation gestalt calculate the restoration 
air content in an average cylinder combining the control unit of not only the first operation gestalt 
but the second operation gestalt. 



[Translation done.] 
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it:iIil360®||A*D<tt/EffiB(c:ffi^C^;(jpAf:i4tcte^-r^yh*!>. JtlEMc iliSWC 

[ 0 0 3 7 ] 
[ 0 0 3 8 1 

^c^:^. ^ (5) Sr^ag-r -5 fcfc o T f:±. ^ (5) SrTIB^ (6) C>ct5t-^?i^LT{> 
[3S 5 ] 



Mci = ^Pmdwnl 1 ^ + wr • A/rf»vw (6) 30 



^ ( 6 ) fC jo T . A t d w n * ffi ^ t m a x ft /Jn M ^ t m i n 4 (7) <D 

B# Tb^ o T . P^&ti(O^T^m^m:ir o i^t^. ^ ( e ) {;i 4o tt ^ n h # ii 

S# t m a X 75^ lb ft /i^ m ^ t m i n 4 T c^) » PbT . ^ /t ^ # 6 <7) P^ # ^ H (il :fcMt ^ 

^ 1^ (7) ;^ n h # ^ ^ gfc » S) ^ v^ t b > ^ hmm^ <d s^'^(0^n\^ 40 

[ 0 0 3 9 ] 



(10) 
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[ 0 0 4 0 ] 

« a IS L v^ i: # ^ tr $ i^ 5 » S L p 
[ 0 0 4 1 ] 

m^:^^h(Dm^mm^m^o kit. ^m:^^>'^(Dm^^Mti.xm.m-r^o ^s^vnt*. :^ lo 

7":y:7'10 2-rMi. ^vVi5':ftir>-i^4 7;^>^^m^(^)iJ^7>'^^CA;ei5^#^i^t5o 
y y 1 0 :^-r5'>^. 1 0 2^;l*5V^T5t#$^^ybr^^7>';5^^^CA35^tB§H#7^^V 

Iv (n) {;iir:yh^tL^flt;6^^ffi^:h.5o ^^fcs SSW#:7^^Vlv (n) ft. ^*Jn 
t-*5Jt5l|imiS(^P^^#60r,S^#t^ffi^^LT:iot?. <D t ^ \Z m i (D ^& % ^ 6 

i: # Id 0 tc-fe b ^ 6 o 
[ 0 0 4 2 ] 

»C V^ T , -7^ >^ 1 0 4 TMl . n - 1 (c: *> it ^ F4 # >^ ^ ^ V 1 V ( n - 1 ) ffi 

:65itj:ir.^yh^tbr*5«91.oBtgiJnfcjottSM«#:7^i/Vlv (n-i) t^ffldso^cir 

h ^ r V> ^5 «g ^ S *JJ ^ ^ tt ^ o i" ?i . ;^ -t^ :y >^ 1 0 4 "C . SfT (U O 20 
^ ^ M # 6 # b T 19 ^ iHl (7) ^ ^ pji SI # 6 # L T v^ 5 ^ o 

^ «9 . ^ (H] (7) ^ B# Pii ^ # 6 # ^ tb ;t ^ T* fc ^ 5 ¥IJ :^ ^ 6 o x 5/ ^ 1 
0 4 :^ V^ T . [Hi O B# ^ # 6 P^ # ^ fc i: # -C *i ?t V^ t 2pJ S S « ^ 
.>^7":yy^l05--i:3itfo 
[ 0 0 4 3 ] 

:^ y 1 0 5 X . ^ n *5 tt ^ g| # :7 ^ ^ V 1 v ( n ) fiS 0 "C fo 5 S 

;^i>BI#^^^rv^7iv^<t*lJ^$tbyh:©'g^ (viv (n) = 0) Idf^. ::^7^5'::^i06-'>«.'r 
5/^li2(^*fT$ix-f.*^{'^7)S^TitLae>t)n6o 

[ 0 0 4 4 ] 30 

-^iiif^. :^'ryyi0 6-^tmi^o ^ y 1 0 6 X . m tf ^ m \c 1 llU :^ h ti 
5o m^:^ ^ ^ ma^ Pii^#6:eiSP^#^tbT;6>feco:$:^{'^(0^tTlHim*^Lrjo 

;^irvi^4 0. mr^'npt-^ 1 9. ^XXf^n^{^^^m±'^-^4 l^^^^ffilW^ixSo 
[ 0 0 4 5 ] 

0 8-1 1 ITMl. 6 mm ^ <0^&nm lhJ^tl(Dm±mP mm a X 

*5J:t;«S/J^fitPmmi n^. S:^ffi^»tma xioJ:t;5g/J^^ii^»itmi ntTSSJSSr^tt 
^ o 40 

y zf 1 0 8"e(^. :^ ^y 1 0 7 X \^ fc ^ Z P m 7b^ m.^^^ ^ inX 

^m^'^P^S.ti<om±mPmma x <t 19 ± ^ f)^ W f)^ . ^J: t> ib \^ ft m M. W ^ t) 

p m^^^R^^ 6 m \^ X h m ± X h ^ t^'s m ^ ix . ^ m ^ ^ t) P m 

t^m±X h ^ t ^ tltcm^ (Pm>Pmma x) \Z (D ^ :^ y- y ^ 1 0 9 Tb^mn ^ iv ^ 

o ;^xy:7^10 9-efi. >^'7^:y7'10 7T-Sl#Lfc®^l^F^J£;tlPm;65P5«l=P^BE:;b(?5S 
:k:fiiPmma x t l^T|2ti$n. m«E(^a*S!ln:d5:g::^^flfBtgiJ t ma x t LTlBH^tbSo 
[ 0 0 4 6 ] 

»:v>-e. ;^'r:y:7^iiOTMi. ;^7^y:7^i0 7T*®j#bfcPii^«f^EE;t;Pm:d5mfiEiaiS 



(11) 
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ffij^Pm^^^S/J^-T^fc^^W^^ix/b:©^ (Pm<Pmmi n) |CCO;^>^x:y>^l 1 1:65^ 
ffi;^<^^/J^1fiPmm i n t LTlE1S^?n. a.*E<OB**Jn;65^/J>ffl^SiJ tm i n i: LTlBig 

$ n 5 o 

1 0 0 4 7 ] 

So ^e^t-. ^ i ^fS <0 6 # L T e> <^ Pii§l« «a^c03»^ffi £ T m ;^ 7" y 

[ 0 0 4 8 ] 10 

ipiSl^^o y h>'W#ii7l^^gf£»m t a v e h ^n^o :iC05(z:i?U;^o 5^ h/U#iiil^^i5ft4 

T*|^^L/cffi;65^iaPii^«^rtM*Tma v e t $:h.So ::<?5^i^P^^l^F^^*Tma v e 
[ 0 0 4 9 ] 

^^v^T-. 1 4-ef:±. :^'r y^ \ 0 9(;l^oV^TM#f$t^/cP^ml^p^I£;^3<7)ft± 20 

ftPmma x^^-r^yV^l 1 llCi3V^TSfT^ttfce5^1^rtjE;^?C0ft/jNfilPmmi n 

Sr®^Uyt:ffi;65®mWrt£E:;^^TfiAPmdwni!$tbS (APmdwn = Pmmax- 
P m m i n ) „ y ^ \ 1 5 TM^ . >^ :y :/ 1 1 1 ^3 V^ T M Sf ^ tl. S /J^ ^ifi 0# ^'J t m 

i n75^e>. >^-r:y>^l 0 9tc:joVNTMfT$ttfc*:^ffl^giJtma x^M^L.fcffi;dSA t d 
wn^^t^-S (Atdwn=tmin — tmax)o 

[ 0 0 5 0 ] 

;^7^:y:/-116-T?J^. ::^7^:y:7'113^115-Tr^ttib^cmtave. Tmave. A 
Pmdwn^dil/Atdwn^^ (6) (;iftA-rSr^{Cj:<9lli^fS<7)^;a^5-^(Dfti 

^{c:yir5/h$^X. PmmaxOffi:;dS^{C, Pmmin(Offi;05ootc, i»^ffl:£mt*5iTJ^ 30 

[ 0 0 5 1 ] 

^r5T-. *:*Jfi?i^fiit^rt^^l3<7)$Jj»igeT*fl. _t5B L J: 5 t;i L Tit ^ $ ix^^ i 

ft ffl] -r S r ^ T t S o 
[ 0 0 5 2 ] 40 
r <^ i 5 k«Si|£fPt#tS^ j*:* Lfc^-&. ^ ISj ffl J- * « ?5 ^ fi o ^ 

fca^). ^Jlsfco^^tc<ti9|g^bTfc'';^h^^3^ifLTtf5^^l-^^-rSj^^^^wu^ 
[ 0 0 5 3 ] 

J: ^9 . mtSrBlT^j^::i^^wU:¥Sr^-(:i^^;tSo. r (7) « ^ S: . HI 5 Sr # P. L T ^ 1 ^ t 50 



(12) 
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I 0 0 5 4 1 
[ 0 0 5 5 ] 

<^ ^ ^ it ^ ^ W T- ti - -t* o o . ^ ^ i§ fci *3 ^^ T ^ i" .5 ^ ^ ^ ^ww =3? ^ ^ 

[ 0 0 5 6 ] 20 

/.C*3. JLfEBftB^-C(i»l^fgSi:m2miiSO^^fiS||cL-ClftPJLfc:6S^ ^T(OM1§i»cov^ 

^-x <D%mmx^ iifsp'^^m^^-SLt^d^ ti^^^i^m (o ^^mmt^ s nA'Xf^m t ^ti. ^<r> 

[ 0 0 5 7 ] 

lll6 5r#KUT. «•^fti^i*^^t5«Si|s|■lSS^fi*JJ;tJ5^ikNF«8Sr^^65«^^^|ll^^:ov^r 
1 2 l^;l*5^^T. «lilHlte®:*JJ; tJ5« H ^ ^ H ig « I- « vn T E c u 3 Wz. ^ 

oT3Sm$tbfcg^3g«Stt:AF t*5<t05iS;^i^^WCA i n j t.jJSBt^^tbS. ifeVNt? 30 

. 7.7- V \ 1 2 x\t^ lll4J;i^Lfcaif^(::j;oTaEa$tbfc^i^1«i-^(^fgrt^E«^«t 

aMc i ;a5^#^tL.5o 
[ 0 0 5 8 ] 

1 2 2{c:*5V^TSt#^;h.fc||iMfg--(;3tt)rt3E«3gm»Mc i ;^7"s/7"i 2 ltc:*^^^ 
r^#$ixfcB«?5;^il:AF tfl^^Lfcffiiis^i^fgilciottS^ilsj-iftltSTAU i i:^ 40 

(TAUi=Mc i/AFt) . «Si|s^lS^^^^^«r<^?^i|s|■l«^^aTAUi<o«s^s^;is^i 

Wj^!iS|-'«lt# 1 1 P}®^ $ ix5o 
[ 0 0 5 9 ] 

^^T^S'-T'l 2 5T-J±, 7.=T y Zf \ 2 2-e^#$tV7t^imfS'--Wfgirt**^^4Mc i 
*>e3. 7=f yZf \ 2 3-C-lt#$ttfca/hffiil^^*3g^*Mcm i nSri^llLfcM^^fJlS^ 
V>T. mi^fSil'*5ft5,'^iA:B#»<oa^^fl;ACAinj i/i5^(±l^tl6„ ''jriJ. iS^4A 
CAinj it±feM5^i:<^li#{*^i6ll^6<)(c:4;t(*ft-^i;:J:oT^tti$n, ECUS 
1COROM34(:1{*:#$^XT*5I9. ;^7^y7°12 5tC*5»t-531^lftACAinj icDgCffl 



(13) JP 2005-90437 A 2005.4.7 

ACAinj i h W t ^ fl ^ o J^:V^t?. ;^x-;/7*126T-tt. @^,^^^^CAinj t 
y zf 1 2 s-ClSm^nytm i i6 ifi, ^ m M <r> m M ^ c A i n j i 

:^ jntzi^i^ . Hi^fSro^^^fflCAinj it^ixS (CAinj i=CAinj t + 
ACAinj i) . :l<OX 0\^i^Xnm^fl1tm i ^n<0^;kmmCA i n j i iC , i 

1 0 0 6 0 ] 

[ 0 0 6 1 ] 10 

®Sl«l^lE;^;-fe:^1^4 0(Z)^l;^^c^^/-rx*s■^*j^-c^^5fc«>, ®^t=rtJE 

-ifeb/i</£ofc!9LTL*5. 
1 0 0 6 2 ] 

l:Sr^»©■^^^'^'/^^cl:oT¥*&^l2Lfc^i (£AT. r5p*&fSirt3te«^^SMc i a v ej 
Ki^^rtlI;^-fe>'f-4 0<^ai;'3lc:y-i'X/i5-^*ixS - tm\^ ^ t) A Pmdwnrofiglci&li 
[ 0 0 6 3 ] 

^ HI o rift -r -5 o Tie *5 X ^ ^^jy 7°i 4 i~;^x5'7'i 5 5*3J;t;5:^7^y7*i 5 9 ft 
. ^ fh ^ fim 4 <D y zf 1 0 i~;^7's/7'i 1 5*3J;tJJ;^7's/7'i i Ttiwi^-e&Sfc 

[ 0 0 6 4 ] 30 

14. *f--1'^'/Ucyc{c:*3ltSf8jrt^*2^M»Mc i (eye) ±iSit (5) icio 

TI1I40;^7"5':;^1 1 6 b^m\^9-m^h.^o 

[ 0 0 6 5 ] 

<^>CV^-C. :;^xs/7°15 8-e(4. TiB^ (.7) (C^Uf::J;p{c:. :*:i^-fi'>'PcycA^e>B'f 
;£IS:Navefc*»t*(l<0-!^'f^'>'V' (cyc-Nave) *»?>*1^-l'^'>'Pcyc*-erofgirt5t 
*^^S:Mc i <^-g-?h$r0f^lS:Na v e-eiSfeStLfeffl:*Si|Z*^eiF*95E«^f5iS:Mc i a v e 

[ Hgt 6 1 40 



eye 

Mciave= ^^Mci(k)J Nave ••■ (7) 

k=cyc~Nave 



(14) 
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[ 0 0 6 6 ] 

fits. m — mmMmTit. m 4 \zm \^tLi^if<Dnt> <o icm 7 iCTik i^tLmi^fA^motx.. ^ 

[ 0 0 6 7 1 

ft p,^^-cv^/j:^^„ iUT, lil8^#flaLT. ^ =. Slife?!^ ffi *i ft 5 IS ^ « 25 w «l ;t * 
fetcov>Tlft?^i-5o 'tis, Ill8»il2i3tl^«|/iia-efcs. 

[ 0 0 6 8 ] 

ti^-C, PAttWl^JE:^:d5|^TU-Cl/^5i:^C0>^D-;/ hyw^jlig^ ^ iSS fi m t d w n 
ir-rsi:, ±IS^ (6) «T13^ (8) WJ:5l-^-rri:*ST?f:5„ 
[3* 7 1 



/{ Vm ) 



JI ;^ <7? ± # S A P m u p i: U T ^ ft -r 5 , 
[ 0 0 6 9 ] 

riDj; 9Jc:A Pmu p Sr3eftLfe»-g-. K^frtEE;^J*SJb#bTV^5»IB]t1'|r*^^t5m 

(Dtitf). ±tes; (3) fiTfs^ (9) <7) X 0 ic^m-r ^ ^ t t^x ^ . m.n.<sp^mt) ±^ 
^m-r ^ :l t tix t ^ . 

1^8 ] 



APmup = 1 •••(9) 



APmup = ^^-^^ mtiq} ■ Alup •••(10) 
Vm 



r«7?S; (10) Sr^J^LTSC (8) {cf^A-TSi:, TIS^ (11) «:#-5rt*s-e#5 
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ISjt 9 ] 



Aid = mtup ' Atup ^. mtdwn » Atdwn •••<11) 

APmup 



[ 0 0 7 0 1 

-t^Ooh. ^ (11) \z Xi^H. l|i^tSi<^®^#64rBI#-r-5ri:(cJ;5Pji^«=[^ffi 10 

1Si'-<^1S5rt ^«^^ftM c i *® Wi- 5 r :6S-t?# 5„ 
[ 0 0 7 1 1 

*jfe-CAPmup*3itJ«A t u p ^mta ' i> ^ t K X '0 . jilM'fc-SiaftirV-l^^ffl 

[ 0 0 7 2 ] 

T»A Pmdwn t . W, i ^M<OPii^#6 ^rj^#i-^mIt'*=Mt5Pii^fF^BE;b<^±#4A 
Pmu p i:lca-^v^-CS i ^tS'^cofSf^^Em^^SMc i Sr^ta-T'S^li: t U T ^5 ;ds ^ 
^ i SvfSiro®^#6 SrBg#-r-5Sfn;:*5tt-5tS^'S'f*9BE;^W-h#SA Pmu p^fttJO 

mi^f§iro®^#6SrW#Lfe«^::*3^t•5®^lSrtffi:fJ<^Jl#aAPmup{cS-3^^T» 
a L T t J; 
[ 0 0 7 3 ] 

5 o 

[ 0 0 7 4 ] 

[ 0 0 7 5 ] 40 

•^fl?cr)ffiStbT3S«)5ri:ibTV^?)„ -f-^jrt?*,, ^ — mmm m ^ if ^ ^ (6) ^^ro 

mt'Atdwn<?3'f^ip"9t;imt- (Atdwn + Atioc)/2^fflV^S„rrt?x 

A t dwn«±aiUfcJ;5(-®^l=rt£E;^:as*:*:ffiPmma xSri:5t#(Da:*:«£^*Jt 
ma xirft/jN'litPmm i nS^i:■5i:§^(75*/J^■(fiB#^Jtm i nteOPeTco^WT-fe'^ (A t d 
wn=tmin-tmax) . Atioc«^i^1^«0K^#6*5 5§#-^5^§iJ 
M) t i oi:®^#6ASH#i-5^*J (W#^«8) t i ciroraro^PBl. -r t> 1b i ^ 
1§<^P<S^#6*SM#bTV>5BfP^-Cfe5 (At i oc = t i c-t i o) „ L-fc^SoT. 

-hlB^ (6) »1. Tia^ (12) ©J;5t-S#lfex.TfflV^bttS. 50 
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[St 1 0 ] 



Met = APmdwn/ \ + mt — (12) 

/ \ Vm J 2 



[ 0 0 7 6 ] 

^(12) -C♦^*^ APmdwnS:^tpJl*S|aiG (b) '^'COEfl|A^cta3■r•S;tf;=^fl:S■ 10 
mtSr^tf«*S|2I10 (b) ^<^SSfBle:tiai-5;*/;^fl:Sr«-t-yi:fe. 
colgirtJEJWS^SMcittiaiO (b) l3©5®Ai:®iBlBi:S:JPSCLfcffii/iot:v^5 

„ HlO (a) i)^ ^ ^ ^ X 0 \' . mfkB\c^m-r^0^^&.^(D:if:^-a:i:-^mt\^XJi^ 

i ffi;^-^t?ii:-C, ?E«^^aSr3R i?) -5 J; 5 U T J: V^. 20 

[ 0 0 7 7 ] 

[ 0 0 7 8 ] 

±oT*fc, l2llltc::^Lfc<t5^j:-^y7'^rLTECU3 1(^ROM34{C«:#:-r^>. ^ 
[ 0 0 7 9 ] 

tC^T*. 6 ^iitili-5.fc«>l;i«^5i:(DPii^# 6 (D^^«J##« ( HI fl?: * -f ) 

^{fbttfc^-^, ■5r^is#^«rt{rfflv^^tvs;^^y :/i/o^Mt;^(-J: •? ecu 3 1 ji-b 

«{c:ffl^^e>^^5*A(D«^^^-i"9 ECU3 l/5^e5^^»#a«'^^B^^^^5g^^^ffi^l 
S^bSriB V^T 9 . 50 



(17) 
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[ 0 0 8 0 ] 

[ 0 0 8 1 ] 

[ 0 0 8 2 ] 

[ 0 0 8 3 ] 
[ 0 0 8 4 ] 

HJfi ?g «i b ^ H H iS ?i5 ^ *5 T li . -9- T' ti5 ^ « ^ m » CO iE ») /J5 /h $ V^ A- :S 

tc i-iii t A. t' fj: m ^ . -t fiioibmmmm^mii^m'^vtmici^ ^m-tizcD nta^titc 

iBi^mt^mmvimi^$} -r^ X /i-mx^p^^m^^&ai'^^^ < mwi i^x i^t. o m-^i-i-t 
. ±mm-mmmmA^ hm^mmm.mi^is i-f :5 J;: X mm ^ tii> ^fH^ 

^ « ^ ^ &/S: ^ SlJ 5 iti il 3iS fo S o 
[ 0 0 8 5 ] 

■So cirofSirt^iK^^S^-r-'i'MlOTli. ^a5-f-.5J;5^-^^>;i^1'^5'/^wf8i^^^«^^ft 
[ 0 0 8 6 1 

tii-r S o 

[ 0 0 8 7 ] 



(t8) 
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[IK 1 1 ] 



^Mci/Ncyl 



/iJo. ^ (13) {CjaV^TNcy l(*5tfS^T$>-5„ Sfc, 2Mc i 1 ^ JU ^ iz H 10 

rtste^^g^SMc i s^l■t^'f^'-'^{;ll:o-c•g•t^-Lfc't)<o-efo5o 

[ 0 0 8 8 ] 

^^»Mc' »cM i «Cl8i tc*fi- 5 MiEfi^l^: u iSrSgLUfeffi^s, ^ i ^ iSi © * o 1Wi ^ 
«^^fl:Mc i' t^^xS (Mc i* = 7j i - Mc' ) o rtt»;ij;i9. #^fSFB^otSrt^ 

^mit^ nm^mtticmn-t^^ t^ti^-r^^ tits. miE^wcj, ^Bg5i«E<^^*s 

[ 0 0 8 9 ] 

#,.^0¥^fflSfcJ*JP«¥^«i:L.TfcJ;V^o mK.ii. miE^^^(OUamW-i^m V i ave 
ttriE^ (14) {wJ:<3T^til$tL5o 
[IK 1 2 ] 



^//ave = S'Tji^„^ + (1 - -s) • - (14) 



7, i(„.,)ttmJ-i^-Y:i'/i'l:::*5V>T3S^ (13) (;i.ioT»i±i$n-:fc«E»^r*fcS„ 
sJiSpm¥*9»*«^T-fof9, 0Ss<i$r^fci-i^*^ft.P3ti'fcfii-C-fc«>o r«J;5(-M 

I 0 0 9 0 1 

;^7^3'ri 6 2~i 6 5*s|I^T$^^>iv^. :^|iafe*tffi*s^«4^fi8H::fc5i:ipj;£$ttyt« 
■g-tcfi, ;5!.-r.y7'16 2'^txitf. 
[ 0 0 9 1 ] 

;^-rs'7'l 6 2~1 6 ffijE-eF.^Tj i a v eOStrd^trtJttSo ^^T^S'T'l 6 2 

V^f. ;:^-r-s/7°16 3(C*5V">T, ;^7^>7'16 2-eSLtiJ$3xfcSi^1S^OfSrt^*2gM 50 



(19) 
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m ^ ti^ o ik^^x . :^'r-yyi6 4X{-is ;^7'-;/7•l6 2^;lt5V^T^iS^^^fcMc i t ^ 
;«.7"y7*l6 3tw*3V>x^m$ixfc2Mc i{cS-:5v^T, ^ (13) {c:j;omiftf§i<c*t 
-t-5MjE#«fc7i i*s|Sffl$ix5, :^-r3/7'l6 5-^tt, -^lHli3±t;5^(Hl;^-r5'7'l6 4»c 
*5V^T^t^$^^;tMiE«^7, ilcS^JvT, (14) l::J:'9mi$ClSft-*ti--5ffljE^» 
(Om&W-i^m. T) i a V e;6S^ai$n-5. 
[ 0 0 9 2 ] 

^fev^-e, :^xy7'i6 6-x?tt, f8irt3fiS«^^fi*7'/UMi0JcJ;oT»W$nfc#*(7? 
¥*S1Sif*3**^^aMc' ;6S^#^tL5„ ■^:LT^■7'-y7'l6 7^;l*3V^-t, xf^-y'T'ie 
5^;lJol/^T^m$ixfcME#^»iPfiipl!D^ii:„ i a v e (;i#5t5(D¥*&tai*i?E*2S^aMc 10 
■ Sr^Jf-rsrfc-e, lli^iaiW4^5fe<7?1Srtite«^^54Mc i' (Mc i' = v 

iave - Mc' ) „ 
[ 0 0 9 3 ] 

ISiF^^JK^^fi^r'/l'Ml Ol'J;oT^aj$*xyt^5l5(752p*SjtS(^?E«^^ftMc' (C^St 

*^^SMc i' (Mci' =Mc' + AMi) . 

[ 0 0 9 4 ] 20 

»AM i it.z(Dmm^mm-r^tzi^<Dh(ox'h^fzi^. :l <Dmm t mm^^m. t ^j: i> x o k 

^if> ^ -^^Wit)^ h ?> . UfciSSoT. f^ffl^VLtM:E;«f:^fiAMit<0M#ttliI13tc:^L 
[ 0 0 9 5 ] 

o. p. q <n X 0 icfs.ism'^m^^^t^iz J: ^ X m-^ ^j:m^ t yj: ^ o 

[ 0 0 9 6 ] 

•^.y:/i:UTECU3 1(DROM3 4(C«:#-r5„ ^LT, 

AMi^aLffli-^o 0'J;t«> ^^ffl^^/65vLlt?fc-5i:#i::aLjll$nfcMjE;«f>^fi>45AMi 40 
[ 0 0 9 7 ] 

-Cigi 3tc;^Lfe-7 2/r*^p>*|]E;Jf;^4AMi;6S|Sm$n, S^i^^Bi'-<^!|f*(OfSirt^« 

^ tl ?> (Mc i' =Mc' + AMi) , 
[ 0 0 9 8 ] 



(20) 
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I 0 0 9 9 1 

m^7^^UM10(cJ;oT^aiStL'5^J?&tgrt?EJK^^S«rMc* . W- ^ ^ ^ M 

*m c ' t i-5o 
I 0 1 0 0 J 

^ ti^?^%(OEEjj) Pat. ±^mm±:^^43\zx^x^ai^^t:ip^mmmmm(D±m. 

3 j;i^3VNr^M ^ttfdP^mSW 1 3 rtoEE;^ i^A^'SP^I^tt) p m t iii^ A :h ^ ti . :iti(b 

[0101] 

Pii ^ I' ^ /l^ M 1 3 d . a h ^ M 1 2 :Jo T ^ tU $ tt a :y b # ii 

S (tJlT. r >pj^tSiF^P^ A^^gEfim c ' J i:«^i-o IS ^ii A ^ ^ SfS S m c * 

^7)^^{;^ov^Tf^. ®^#^7=^/^'M14^C:fcV^Tpa^-r^) i;65A;':$tL. ^ X t) ^ 

3 :}s XXf^ — ^ :^ ^ 1 4rt^c#^Ei-5eSm;^:^<oj[E;^; iSXT. rpii^@rtJBE;^P 
mj ir ® 1 3 *D J: V? ^ 1 4 S K^;^f j:^ (DMS (WT 

. rp^^l^rtMSTmj ir^^^-T) t^^Stti^tL^o p^m«^7^/i-Mi 3^^^3v^T^ta$n 

fcP^mffl^EE;^?PmtPiiMl^rt?aa[Tmfi*(CBS*#^7^yUMl 4--A;^^ix. ^e>{CPii 

\^&t} P m{t :^ t3 y h^w^-r/UMl 2(c:tA;^$ti.'5o 
[0102] 

Pii^#-=e7^/U'M14t;i(^. ^An^ JUM 1 3 \Zi6\^^X Wtti ^ tltz^R^^ [^J^tl P in:}6 

XU^R^mi^umTTn(Di&\z±n,umT a ;65A;':;^ix. r^^t.A:;':;^i^y^c^/^'^^-^<^ 
ffi^Si£E'r5Pii^#^7'/^M14<7)^T^/l-^{c:ftAi-5rtT*. ^iSfSF^KA^SliJEam 

^i^^a^n^o ^ttl^n/c^i^f§5f^?gA?^^«D»mc' it. ^i^m^^m^^MM c 

' \z^^^ti^ ^ (Dwmm\h^m^m.mM c ' {zm-^\^^xmnmM^^^h(Dm'^mmm:d^ 

^^^tl^o ^tz. mm,Wey'/UM13\C^\,^xmm^tLtiW'^m^^RA^%i&ATnc' 
PiJ m t= ^ Ml 3 ^tb^o 
[0103] 

mi 4. h ^ ^ X {z ^ Wi^^m^nM^ef/i^Ml 0TM^^^^x>'V'Jc::^oV^T^m^ 
tLfc/^•^p<-^^0ffi;6SSlJO^7^/u-^OA;^^fli:^:bT*^Jffl$^^5<^T'. 

/\^MlO±i^Xit. mm\z At} ^ ti ^ m.it :^ ^ V h /\-m A e t. ;k:^JEPa. ioJ:!;^;^^^ 

»M c ' ;d5^ttl ^n^o 
[01041 

25c jc:, fSp^^iiE^^fi^x/HV[l OO^^-t^/u-MI l-^Ml 4(cov^TUie|■t-'5o 
[01051 

S;^®JW:^n-/h/u^7=^/UMl 1(1. ft?tir>'l^4 6tc:J:«9<^UJ$ttfcTi^iryi-^^/u 
Slf^S:AccptcS-:5V>T. ^^o:^ntyh>^u#18;5S3f^^^^AT^Jcaia-r5;^city 
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fzT ^ ± /\^^ y /u^i^m A c c p t s m 1 5 i^tzr ^ ± /i-^y /i^^i^MA c c p t n 

bti^ o w<o±5J::UT3^^P>*xfcxn-/h/vB§a[0 t n :^ ^ h ji^ ^ /^M l 2 ^ t tS^ 

htl ^ o —io. ^ <D u y Y jvm^ Q t ^m'^^m ^ T 6 4msec) ^h'tt® 

[0106] 

r O J: 5 . g ^ n h /u- M S 0 r , ^ ^ ^ ^ fi^ ^ 0# PpI A T giT tTD ^ ^ *5 
5r^5^^r/^-<^5^/u^f^4Accp^cj^;;CT*^$^^:rc;^oy^/ugBS0ti:^bv^e iH 10 
:?^ny hyvBISe ^^^:S<5v^T:?^D-/ Y jv^ \ B ^> mWl $ ti ^ tz ii> . :^ xj ^ h /\^mm 6 t 

V> 5 o iS? I- M 9 t . :^ n :y h S 0 t , H >^ n h /V # i 8 @ff ^ ^ PbI A T ^ 
PJ ^ -r 5 :^ D h S ^ 7i o T V^ 5 o 
[0107] 

m :^ ^ y h/v-^T^/UMl 1 t)^ (b iti tl ^ tz ^ :^ ^ y h m 0t75^?^.Tfa^(15 
) ^3:S-^V^TP^cl-/^/^#iiia^M^5fi»Int;6S^m$n^o rrr*. ^ (is) tC:}bMt5 

^ L fc J: 9 >^ 7^)^ ?3 * :S ^ o * . A t n h /v # ct) 5fi p ^ B ^ §r ^ b . >^ n 20 

:y h /V- * 0 t <Z) 11 ^ . 11 1 7 ;^ L <t 5 :y >^ ^ * ^ o ^J: ^ . :z iv h ^ 

S^^M^Bit/gHPfgrffi^AtSrSirdbfciM •At^:^n>;/hywr^Se t;0^ib — oO-^^;/, 

1 moia t) (^m) <^SfiMlmol-C^S[ LfcfliT^feSo 

[iS 1 3 ] 



[0108] 

^fc. O (Pm/Pa) (ITIB^ (16) \Z tt^ fz ^ . ^ (D ^ (16) {;i^bMt 

^/cfiit^lt { — -^mt-t^) T-fe^o r<^II^O (Pm/Pa) til2l8lC;^L/c:J:9 

3 4lC{ft#L. (16) ^m\^^X%^W't^(OX\t.t^<'^yZff)^h<^ (Pm/P 



(22) 
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/c-l 
2k 



2{k + l) 



Pm 
Pa 



1- 



Pm 
Pa 



Pm ^ 1 

Pa ~ *r + 1 

Pm 1 

> 

Pa /r + 1 



(16) 



[0109] 

:i tl h :^ o h ^ f f^'M 1 2 <D ^ (15) JBit/SC (16) li. x^ns/b/U#18± 
iSSro5ci!|s:rojE;^3S::*:55JEPa. :^ u y h \ B ±.m<o%W-(r>m.$L^ -X^mSlT & . 

y h # 1 8 » ^ 7^ *f L T . « » i*: # lU . /u =¥ # MiJ *5 i l>* il » » # MiJ ffl 

[01 1 0 1 

• , *5J;t)f:fe^?a*Ta*^p,. TI2S: (17) iait/SC (is) icS^v^•CPiit^l=rtm;'3 

Pmfc-J;t)5®^l=l*9taaSTm*SSLtt!^ctT..5o ^Ct5, ^ (17) *3j;t)«^ (18) {C*3»t5 
Vm(4Xnyh7W#18/4^bP^^#6*-C<73K^Wl3^<OSlJ^^ (£AT, rP^St^flFB^J-J 

i:?4;-r)<7)Sffiii^Lv^:fe^-efc5o 

[|^ 1 5 ] 



d(Pm\ R , 



(17) 



dPm Ra , ^ „ ^ 

— — = K {mt ' Ta - mc-Tm) 

dt Vm 



(18) 



[0111] 

»> i-'5eip*>^*&fgii*9KA^t(«ftfim c ' t w^tc^ Lv^fc^^), «&«#IiJ i •? Tffi* 

(19) r<D^(19) *Jj;U!^flC05i^Jg3^gS;(Pm-Vm = M-R-Tm 

) iU. ^ (1 7) 
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[l» 1 6 ] 



= mt-mc —(19) 

dt 



[0112] 

m^v&^'SP^m.mT m t -t ^ t . /^^^^wmK ^ 'o rfs.^ (20) i^n^ti. r^s; 

(2 0) t5J;tJ«±IB^tt:ro<^lS*m^J;9. (18) ;d5#p5tu5. 
[SCI 7 ] 



20 



30 



^^^L^^:l!^=CpmtTa-Cp-mc'Tm -(20) 



I 0 1 1 3 1 

iS^#^7^/'V'M14-t?tt. PJi^l=rtEE;'3Pm. PJ5^«=rtMaETm. *5<tt>*:;lc^?l.STa 
(2 1) |c:*jft2)a. b{4.«BllHlte^Ne;ir>e,. Se)tc:K^#6<^fiC*B^ (>'-<>'W:/i?-r 

1 8 ] 



mc = —iaPm-b) —(21) 



[01 14] 

Mc' J4, ®m«=rtjBE;t;Pm(c:jt0iJ-r-5i:iai£ti-5c:tAs-etSo 

[0115] 

(7)^1g(0®^^fa^^:aoT (^ifii-5J;5l-*llifi?i^^-C'tt^7>'i?^i8 0° :») 

bfe t JPJ^ISI^® A^^SSSm c ' (UA TT*PaE-r 5) i: i" * i: . ¥ ^ fS ^ « ^. 

Wl*3JE>3 P mtcjt^tf-r?, t %;t bixSo r o r i: > 3ggS *5 <t Otg ^lU »c * ^ v^ T . ± 50 



(24) 
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lE^ (2 1) i5#btt5„ ^ (2 1) (Cfctt'SffiaJlit^J^IS-Cfc'J. ItlilUte^ 

^''-C^*»^^fcTa/Tm75S«^^$i^TV''5o 
I 0 1 1 6 1 

foSli-g-JJiov^TlftW-t-So 'fcfc. m 2 2 7 ^ y h <D M.m fkBt^ KE W *s 

55^ftfl 3;4^ib±^f^(D^j^^lcgftAi-5 5S^;>?/;^»Sft*(il2l2 2(c^^T-:T^Lfcil!P 
•rfeS. Sfc. 0iJ;t«l#^fS'^<7)¥*&f§5F*33E*§g^»Mc' {ill2 2{;i^l8-C^L7tii 
•3 f fe 5 , 
[0117] 

rttlw^fLT. ll«8t?*bfc®»W*^b^T<^ft1^ro^j^^tciJfDA-r^?5$t;!/xwftSr 

r <7> — ;*C«i6S-eS^ L fc¥i^fail*3® A^^«Efim c ' . 4 M fS » til »4 ^ 7 ^' -Y 7 
hd^i 8 0° (■*-''£*:>*>> 4><; ho-^^|^^^l^^c*3V^Tl■i^-l'^/^=^^;I;^9V^'v■^ 
7h;dS[Elte-r5^«7 2 0'' 4r^m»t?fiJofc^lS) m ffi i" -5 <^ I- S ^ Pel A T 180° 
Sr^lfLfc't,<D?5Sip*SfSi(*35E*3g^SMc" i:/i?.„ LfciJSoT, Si # * x M 1 4 T? 
ffai$ttfc3Fj$ifgif^P5A§^^ii!Eamc' {C A T 180° ^^^-TSitr-, ^*&1Sirt?E*^ 
M»Mc' ;!i5^ta$tL-5 (Mc' =mc' • ATI 80° ) „ J;i9l¥;isfflt-«. W-^^^%m 

!RA^^gE»mc' »c A T 180° **SLLyS:tW*S¥^iart*«^«t4Mc' t^tvSc 
[0118] 

f^JE^^^^i^T^/l-lOSrfflV^T. ±|S^ (15) , ^ (17) , (18) , tJitJf^ 
(2 1) ^fi?<ii:(CJ;'9^^tt5o i£0:^-g-. ECUSl-C^O-S-TS^tol:!. CDjxf^c;? 

^^fiSticfb-r^itJ^^^dSfc^, ^Sijt, ff^FB^RSA t^fflV>T;* (1 5) . ^ (1 7) . ^ 
(18) . fcitZ-S (21) SrfittfeYki-Si:. ^tb^HTfEiC (22) , (23) , 
(24) , ioiO?^ (25) ;OS#<b*xSo ®mt=f^faSTm (t + At) {4, S: ( 

2 3) (2 4) l'J;oT^^x^*tt»tii$ti.fcPin/Tm (t + At) fcit/Pm 

(t + At) t^h . (26) (;:J;oT#tti$tT.-5o 



(25) 
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(f + A() = ^(0 + AI • ^ . („,(,) _ „(,(,)) ... (23) 

Jfn I fit vnt 

+ A/) = /'/w(r) + A/ • *: • ^ • • Ta - /wc(0 • Tm{t)) • • • (24) 

Vm 

mc{t) = • (a ■ Pm(0 - 6) - (25) 

7m(0 

Tm{t + A/) = - — - • (26) 



10 



20 



[0119] 

(7)^ (2 2) r*^ffl^ixfc^Satt;i;*ott5^n:yh/l^#iligi^^?5ft4mt (t) PJi^ 
#^7^/1- Ml 4(^^ (2 5) -C»ta$fbyh^^JttC:Joft^^i^^F^PAA^^«KSmc' ( 
t) t . m^"^ f JVM 1 3 (D^ (2 3) :JoJ:t>*^ (24) \- iX K ^ ^ . 
^Jt + At{;i^tt^P^^WF^i£:;^Pm (t + At) *5J:t/PSmi^f^?aSTm (t + At) 
t^^m^th^o ^feV^T% »tH^tl.fcPm (t+At) ^SctU^Tm (t + At) 

fyU^'7^/UM1243<ttJ^Pii^#^7'/VM14(7:)^ (2 2) *5ctrJ^^ (2 5) (Cf^A^tt. 30 

:itb{rj:i9Bt^Jt + A ttcfctt^:^^-/ hyw#ai§?g^»fL»mt (t + A t) -io ^xi^^ 
1SrtPiiA^^»ft»mc' (t-fAt) ;6s»aj^n^o>bT. ;icoJ:5?i^f+#^i^»:)3g-r 

co^^j t tt -5 ASg^gfiim c ' ^s^ffl^H. St ffl $ tt ^piS 1S A ^ M 

?5S»mc' I- ±|S^P«1 A T 180' ^ S ^ i" 5 r ^ . ffi S B# t *5 ^ ¥^ ?5 J* 

^^*Mc* ^5»m$ix^o 9tm^:^ ^ Y jvmm Q t;)5^|^(7);^n:yh/w#18 

[01201 

^£Ei:^UV^ ( P m (O) =Pa) i: $ . PJi M * ft * T m M M S ^ ^ U ( T m 
(0) = Ta) t^ttT. «-^T'>'l-Mll'-M13JC*5lt5ff-»;&5rai6$tt5o 
[ 0 1 2 1 1 

/^fc^. ±fafiirt^«^»a^x/i-MlOTMi. ;'c^MffiTa*5J:Ut:^^ffiPa;()s-^T- 
^ b T v> 5 >&5 ^ Nf^JJc J;oT^<ti-Sffi^ t <t < . X. . :^ ^ S «r It tti ^ 
/c:^6O:A:^M|^^rvi^{cJ:or0#^Jt^c:^oV^T^tW^Hfc^itSr::^ (t) . 

JES:J^aji-^/ca6<D:^^iE-fe>'i^(;:iJ:oT^giJ t \^^\^^ X ^ fx^t^m,^ & P a ( t 
) i:UT±fB^ (22) . ^ (24) . ^ X.-a^ (25) {;i{tAi-SJ:9l-UTt<tVNo 

[0122] 
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■7 



-9 («^Jx.«. 

5 «c M 5^ ^ tt r 

2 3 ] 

1 f =7^^ if ^ 
] # fvf^ 

6 1 # St fiS tc: *3 

fe 5 „ 

7 ] #^f^ to 

- h T' fe 5 

8 1 11 = 2 

9 ] ^ 185 to Pii 

-ria t? S 



1 



1 0 

1 1 

1 2 



- h t? fe 5 



HI 



P 7 i/ 

H i ^181 



X D h 

MiSc <& ( 
X n y h 



Blc7^$^J^ft$lSJJ5fflv^P>^^5 ^ « 53 S% f 0 T? 5 o 

o 

o 

Pm/Pa) ^^-f-UT-fc^o 
'7^/W(DS2|s:|R;g:S::^-t-|lIT*fo5o 



10 



20 



30 



^ CO tft ] 
12 4] 









5 






6 






7 






8 






1 


1 




1 


3 




1 


4 




1 


5 




1 


8 


••• u -/ h yu# . 


1 


9 


:ai T 7 n p< — 


3 


1 


-ECU (m^sij^ai::!- ^ b ) 


4 


0 




4 


1 





40 



(27) JP 2005-90437 A 2005.4.7 



[Hi] 



im 2 ] 




LS I ^mni'g | -i — |A/D 




ml 



Pm 
Tm 
Vm 



mcl 

. mc2 
_«=> 

- mc3 

- met 



Pm : laSlfl'rtffACPa) 



[IS 3 ] 

(a) 




tmax tmin ^^i/OM<iBg) 



im 4] 



Bt (start) 
I n-n+1 




^113 



107 



Pm.mt.TmgC^ f -^ 
108 

Pni>Pmmax ^"^--J!:*o 




Z mt*- Z mt+mt — 112 
ZTm«-ZTm+Tm 



mtave«- Z mt/m 
Tmave*- Z Tm/m 
I 



APmdwn 

«— Pmmax-Pmiran 



Atdwn 
—tmin- tmax 



115 



II Mciyta 1^ 116 
117 



m^«-0 
Pmmax*-0 



Zmt— 0 
ZTm— 0 
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CA1CA2 



TDC 



[IS 6 ] 



EB6 



(start) 



AFt 

CAinjtiRW 



^ 121 



I 



TAUi 

^Mci/AFt 



X 



Mci-McminJ: U 
ACAimmm 



^125 



CAinjj 

♦-CAinjtt-ACAina 

C end) 



— 126 




(29) 
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Mci 



ZMci/Ncyl 



.165 



II I? iav e— S- n i{n)*( 1 -S) t? Vb- i > |] 



i 




VL1 



VL 



(30) 
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im 1 4] 



au 



4 



E 












^ 








CD 





K 




Accp 



[ IS 1 6 ] 



B16 



ju(et) 






Pm : flftaertff *(kPa) 
Pa : :*:ftff(kPa) 

Ta : :*:SliajS(K) 

At : :^ay h;u#gaDtBTffia(m*) 



0.2 ac 0.6 

Pm/Pa 
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2 0 1 
a 20 



mt. Ta 



J 



Pm 
Tm 
M 



13 



im 2 2 ] 



1922 



mc.Tm 



Pm : iafttt(*3E*(kPa) 
Ta : WLX1-i»SLi^<r)M&(}0 
Tm : Rfttf F*3aA(K) 
M : Rft9l^(0»ftft:S 
R : ftf4:£tt 
Vm : (RflC«S«(m*) 
Cp : ^fBLOy^Zttm 
Cv : ^S<7?3g$it»ft 
K : ttmit 





Pm : flRfi#J:sa©ff*(kPa) 
Tm : ®Sl#±35E©aS(K) 
Ta : i^SLiSSCK) 
mc : eSF*9?&A:tf;^3i£S(g/seo) 
a^ : 



(32) 
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(SDlnt.Cl/ FI r=i— K(##) 

F 0 2 D 41/36 B 
F 0 2 P 5/15 A 

F^— 3G022 AA03 DAOl DA02 FA02 FA04 GAOl GA06 GAG? GA08 GAll 

3G065 AA04 CAll DA04 FA12 GAOl GA05 GAIO GA26 GA27 GA41 

GA46 KA33 

3G084 AA03 BA04 BA05 BA13 BAIT DA04 EA05 EA06 EAll EB12 

EC03 FAOl FA07 FAIO FAll FA33 FA38 

3G301 HA06 JA03 JA18 JA20 LA03 LBOl MAll MA12 MA13 NA08 

NC02 NC07 NCOS ND02 ND06 NEOl NE06 PAOIZ PA07Z PA09Z 

PAIOZ PAUZ PEOIZ PE03Z PF03Z 
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[%tT0l 18^8>^ 24 0 (2006. 8. 24) 



I M # -i- 1 2005-90437 (P2005-90437A) 

^^17^4^7 0(2005.4.7) 

[ tBS#-§- ] ^^Pli 2003-327965 (P2003-327965) 



F O 2 D 
F 0 2 D 
F 0 2 D 
F 0 2 P 
I F I ] 
F 0 2 
0 
0 
0 
0 
0 
0 



F 
F 
F 
F 
F 
F 



D 
D 
D 
D 
D 
D 
P 



45/00 
9/02 

41/36 
5/15 

45/00 
45/00 
45/00 
45/00 

9/02 
41/36 

5/15 



(2006. 01) 
(2006. 01) 
(2006, 01) 
(2006. 01) 



3 
3 
3 
3 



E 

J 

D 
F 
Q 
B 
A 



immB] ^^18^7^ 11 0 (2006.7.11) 
[^j^MlE 1 ] • 

im^m 1 1 

m t . 

. ^ ^ fS *f /Jr. i- ^ PA ^ # ;ds 



i^-r ^ ^ t X \^f^t) go M ih # 

imiEM^m @ ] 0 0 0 8 

[ 0 0 0 8 ] 



0 & T » 



^ n h ii il 
^ # 53 # -r a 

T «l a ^ tt :^ 



(2) 
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-i^mmiE 3 1 
imiEn^^ s^] 0 0 15 

[ ffi IE <73 5& 1 
[0015] 

%'s^mt)i)^myh k fj: ^m^bm.^m t cor^(Dmmi: t dwn t l.. jifsp^^#(OBs#^ 

«tM#^«toP^rofflF^SrA t i octL. r ix P> ffl »c v^ r Jn IS SS 4 -fe :^ iMc J; 

m#®tt. n y h /uiaiS^^fiM t Sr TIB^ (1) Ic S -:5 V ^ T » ffl i" 
Mt=mt' (Atdwn+Atioc)/2 — (1) 

imjEM^m g i& ] 0 0 3 3 

[ 0 0 3 3 ] 

^ (3) ±19. P^^«lfi5 5)---<^«SX^aM.* (mt) /JSgitm^^SfEft (mc i) <0 h 
±^ l'fthitv&n'SP^S.t)*^±^L.. /J> $ It n « ^ 1= 1^ EE ;'3 ^5 1^ T b . ^ b it « ® ft 

h^m^. o&^'sm^A-hff^mm^^mm (mc i) ttm.^^ e <Dmm\^^^\^xmx6^ 
x'ki>(D\zM p&%'§^^7!)^ h o^mA^nmrn. (mt) cxv^ ^ mm t ^j: '0 ^ <o 

^itt^m-^t^X ^ . rrofc*. «5Etti25M«fEa (mc i) taEAS^iffifi (mt) tO:^ 
•? jS-T (ma (a) #fii) o rtL«. ±fE^ (2) (0^mrc>*s/='f*3<0 

imiEn^m i i& ] 0 0 4 2 

[ttjE:^fel 3gM 
[Wioi^?^] 

[ 0 0 4 2 ] 

i^V^X-, ;=^X-y7°104T-»4, ^a"Jn-l{c:4bMt-5^M#77^^Vlv (n-1) Wlli: 
1 iC± h ^ tlX is y) K^mm n \CtS i-f itmmit'y 7 if V l v ( n ) m.*^ O \Z± iy h 
^ fl X ^ X h ^ f)^ m ^ ^ ti ^ o -r^i >^7^-;/7°10 4T'{i. SuHK^^f^ 

*3 T . <^ iHl ^ AS K St # 6 AS ^ # $ Ix i: # -e « i: fij « $ tb # ^ . 

T^.yT'l 0 S'^irittP, 
[ ^i^ffijE 6 J 

imiEnmm s ^ ] o o e 4 

[ttiE:^jfe] ggM 



(3) 



JP 2005-90437 A5 2006.8.24 



[ ar je © ^ ] 

[ 0 0 6 4 ] 

^'tt, «llis<EM^6*»e>w-!^-f^'/^»^«-r*<>^^^'T-fo■5o ^^:v^-e, i^^y^i 57-e 

. 2|s: 1^ ^' -'^ c y c i:: Jt 5 ISi fiM c i ( c y c ) *s , ±|BiC (j6_) J: o 

T m 4 (7) ;^ X y 7" 1 1 6 t 1^ «Hc 3¥ m $ 5 „■ • . - 

[^i^^tE 7 1 

tMjE*t^3SS«] 0 0 7 5 
[ 0 0 7 5 ] 

^?i?(OS«i:LT*«)6ri:i:LTV>S„ "T/it^-^. ^ - ^ ig fl^ ^ (C ±3 tt 5 (6) 4^ «0 
mt • Atdwn<DftiP>9lCmt • (Atdwn + At i oc) / 2SrfflV^5o r^t?. 

ma xi:*/jNffiPmmi n b. ^ b ^ (O Mc 'h M t. m i ntroWcDBtFel-efot) (A t d 
wn=tmin-tmax) , Atioc«||i^fgiCO®^#6;65^#-t-5^^J 
n) t i oi:®5t#6aSH#-r5^»J t i c b (Dm<r>^fSi. -ttiti^m i ^ 

fBroK^#6;iS^#bTV^-5^F^-r'fc-5 (At i oc=t i c-t i o) „ Lfc^iJoT, 
^Eg^J£Ji^^r-tt. ±IES; (6) fi. TIBS: (l 2) ^DJ:9^;l##mx.T.l^^e>^^5o 



[ 8 1 

[ If E a ^ 1 0 0 7 6 

[ ffi jE CO F*9 ^ 1 
[ 0 0 7 6 1 

^(12) -eti. APmdwn?r-g^|>3ld5|3l0 (b) 4"<^ffi«A(C+i^-r^;tf;^»5r 

«fSrt^*3g^SMcifi:llllO (b) 4'roEi»AtffiffiBi^^JP^Lfc1i:t'ioTV^•5 
, laio (a) *»P>5J-!5>5J:5(c:. ffi 81 a_{c: +i ^ -f- S ;<f S lit ^ 1: •& ?^ L r * 

ir^'ioTiot), #^iai©K^#6 LT 5 fc oT-Jifflirt^*^^ 

bm.^'^t>'^X ^ Jp*&faSl^*JK^^aSr*feS± 5 ^-LT'{>J:v^„ 



[i» 1 0 ] 




•••(12) 



